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applications in the field of aeronautical endeavor. By including in its membership 
the leading aeronautical specialists in other countries an interchange of international 
aeronautical scientific progress is made possible. The organizers of the Institute 
believe that the future development of aircraft depends on scientific research and engi- 
neering experiment applied to design. As practically every scientist and engineer 
already belongs to his own special professional society, in order to bring together all 
specialists contributing to aeronautical progress, the dues of the Institute have been 
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Longitudinal Stability of a Bi-Motor Transport Airplane 


CLARENCE L. JoHNson, Lockheed Aircraft Corporation 


(Received January 11, 1935) 


INTRODUCTION 

“T*HE problem of longitudinal stability has received 

considerable study, but in spite of the volumes 
written on the subject, the practical designer is continu- 
ally running into difficulty with it on almost every new 
design. This paper deals with the author’s experiences 
with stability in the development of the new Lockheed 
Electra bi-motor transport airplane. Having carried 
the problem from the wind tunnel through flight tests 
and finally the Department of Commerce tests, he feels 
that, while not adding appreciably to that already known 
about longitudinal stability, a survey of the problems 
encountered might help others avoid them. 

The reasons why an airplane must be stable are more 
or less obvious. A moderately unstable airplane may be 
flown safely if the pilot continually manipulates the 
controls. However, this procedure is very nerve-wrack- 
ing and tiring. Landing or taking-off is dangerous in 
an unstable ship. In flying blind, it is a great comfort 
to know that if the controls are released, the ship will 
continue on its course safely even in rough air. It is 
possible for a ship to be so unstable that the controls 
are insufficient to pull it out of a dangerous attitude in 
flight. 

The use of automatic pilots (such as the Sperry) will 
probably allow slight instability in the cruising range 
in the near future, but stability at high angles of attack 
must be maintained for landing conditions. 

The following discussion deals with longitudinal 
stability, first in theory, then from the point of view 
of wind tunnel and flight tests. A discussion of the 
theory must necessarily be little more than a statement 
of the factors involved, in a paper of this length— 
however, as a review of the facts, the following is 
offered. 


THEORY 


If an airplane is statically stable, aerodynamic forces 
must arise which tend to dampen out any oscillation in 
the flight path. Static stability is generally considered 
apart from dynamic stability, which deals with inertia 
effects in flight. Only static longitudinal stability will 
be considered here. 

In the conventional airplane, unstable moments arise 
from the wing, which must be overcome by stabilizing 
moments from the horizontal tail. The important 
variables which enter inte the stability picture are these: 


(1) Location of the center of gravity of the airplane 
on the mean aerodynamic chord. 


(2) Wing chord. 

(3) Tail area and aspect ratio. 

(4) Distance from C.G. to C.P. of tail. 

(5) Blanketing of airflow over the tail by the wing, 
fuselage or nacelles, 

(6) Propeller slipstream. 


The factors are dealt with quite thoroughly by 
Warner? and Stalker?. Professor Stalker has derived 
an expression showing the effects of variation in the 
main variables on the stability criterion. This criterion 
is taken as the slope of the C,, vs a or the Cm vs Cz 
curve for the complete ship. If dC,,/d C, is negative 
throughout the range of angles of attack, the ship is 
statically stable. Where C,, = 0 the ship is trimmed 
for flight with the given stabilizer and elevator setting. 
Neglecting propeller effects, Professor Stalker’s expres- 


sion is: 


1 Edward P. Warner, Airplane Design, McGraw Hill, 1927. 
2E. A. Stalker, Principles of Flight, The Ronald Press Com- 
pany, 1931. 
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TABLE I 
| = 
h At (Ry-1.3) | dC, | 
Type of Plane 371 Aw Ry < (R, + 13.7) dC; 
| 
Twin motor flying wing — 15,000 | — .0168 . 1665 .80 1.985 | —.0567| .97 
Low wing monoplane, fillets, re- | 
ee | .342 | —.0258 


High wind braced monoplane. No! 


| 
236 | .519| 2.12 | 1.072 | —.091 | .85 
| 


1.330 | —.171 | .744 


181 .450 2.70 


| 


Low wing, 2 place cabin monoplane, 


| —.109 | .937* 
.932 —.061 | .708T 


.300 | —.0149 .125 .630 | 2.84 
High wing tri-motor. NACA cowls, | | 
| + .0280 | .470 | 2.54 | 1.200 | +.015 .48 
Low wing, retractable gear, no fillet, | | | | 
| 332 | —.0218 | .177 | .651| 2.46 | 910 | —.212 | .744 
Low wing, retractable gear, Prestone| | | | 
| —.0227 | 186 | .677 | 2.46 915 | —.167 | .910 
| 
| 
Twin engine, retractable, low wing. | | | 
.330 0 178 .706) 2.73 | 874 | —.130 | .770 
| 
* Vee shaped windshield. : 
+ Flat, rearward sloping windshield. that ships which showed a value of dC,,/dC, of —.15 
or less (as —.20) by the above formula had satisfactory 
dC,,/dC, = — 23 4 a/e — B7hC,/c — (Ay/Ay) longitudinal stability in flight. It is recommended that 
xX (Ry/Ry) (1/e) (Ry — 1.3)/(R,y + 1.37) Je (1) fast transport ships show a value of dC,,/dC, at 
least as low as —.20 in the cruising range. 
where 
a = distance of C.G. back from the leading edge 


of the M.A.C. 
c = length of the M.A.C. 
h = distance of the C.G. above (—) or below 
(+) the chord line of the M.A.C. 
C, = absolute lift coefficient at the angle of attack 
being considered. 


A, = wing area including fuselage—wing inter- 
section. 
A, = horizontal tail area. 
R, = aspect ratio of the wing. 
R, = aspect ratio of the horizontal tail. 
1 = distance from C.G. to C.P. of the horizontal 
tail. 
e = the tail efficiency factor. It corrects for the 


frictional loss of the airstream in flowing 
over the fuselage, wings, etc., before it gets 
to the tail. (See Table I.) 


The use of the above equation is chiefly as a guide 
in the study of proper values of tail area, tail length, 
and C.G. positions for a new design, or in the study of 
the wind tunnel moment curve. The equation is appli- 
cable mainly to the high speed or climbing speed range 
of angles of attack. It has been found from practice 


CL 


08 of Attack 


fo 


Fig. l. 


Fig. 1 shows a typical wind tunnel curve of C,, and 
C, plotted against angle of attack for an airplane model 
derived from the usual tests, with fixed elevators, with- 
out propellers and taken at a constant wind speed. 
These conditions by no means simulate flight but up to 
the present time, the wind tunnel still furnishes the 
most direct way of studying stability. The product of 
(dCm/da) X(da/dC;) = dCy,/dC, as referred to 
above. The decreasing slope of the moment curve at 
high angles is caused by a blanketing of the tail by 
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the wings, that is, by a decreasing value of the tail 
efficiency factor “e.” “E” may be considered as that 
part of the energy in the undisturbed airstream which 
is recovered at the horizontal tail. For present day 
aircraft “e” varies from .75 to .95 as determined by 
taking the slope of the wind tunnel moment curves and 
solving Eq. (1) for the value of “e.” In this manner, 
“e” includes a correction for various discrepancies in 
the derivation of the formula but the method is thought 
to have some value. This was the method used in 
deriving Table I. 

The most effective way to increase the slope of the 
moment curve is by: (1) moving the C.G. ahead on the 
M.A.C; (2) increasing the tail area or aspect ratio; (3) 
decreasing the wing chord. Quite often an investi- 
gation of windshield design, nacelle location, strut angles 
and so on leads to surprising increases in stability as 
shown in Table I in one case. 

It can be proven theoretically, and flight tests sub- 
stantiate that dC,,/dC, varies linearly with C.G. 
position. German tests* show that 


[d/d (C.G.)] (dC,,/dC,) = — 01 
if C.G. is expressed in % of M.A.C. 


PROPELLER EFFECTS 


The effects of propeller thrust on longitudinal sta- 
bility are illustrated qualitatively in Figs. 2 and 3. Fig 
2 shows the case of a simple wing and tail combination 
and illustrates the effect of the wing downwash on the 
angle of attack of the tail. With the propellor effect 
included conditions might be represented by Fig. 3. 


Power off condition 
Fig. 2. 


Va 


~ Vp = Vy 
Power on condition. 


Fig. 3. 


*N.A.C.A. Tech. Mem. No. 708. 


Fic. 4. Model of the Electra in the Wind Tunnel. 


Another effect is present if the thrust line does not go 
through the C.G, An increment is added to the moment 
curve equal to the change in thrust with angle of attack 
times the distance from the C.G. to the thrust line. 

The effect of the propeller on stability is so complex 
and varies so much from ship to ship that it is seldom 
calculated. 

Winp TUNNEL Tests oF THE LocCKHEED ELECTRA 

The wind tunnel tests on the Lockheed Electra were 
run at the University of Michigan on a model of 55 
inch span at a Reynold’s Number of 545,000. As the 
model was first tested it had a conventional single verti- 
cal tail, a large wing fillet and the horizontal tail was 
located as shown in Fig. 4. This figure shows the 
motors removed and the nacelles are completely faired. 
The moment curve for the original ship is indicated 
in Fig. 5 as Test 1. All moment curves were com- 
puted for a 33% C.G. position on the mean aero- 
dynamic chord. In the cruising range, dC,,/dC, was 
—.089 but instability and neutral stability are indicated 
for a large range of angles of attack. An investigation 
was immediately started to determine the simplest 
method of correcting the above. It was quite apparent 
that the horizontal tail was being blanketed by the 
engine nacelles especially at the higher angles. The 
engines were removed and the nacelles faired as shown 
in Fig. 4. An immediate improvement in the slope of the 
moment curve was noted. (Test 2, Fig. 5.) An increase 
of slope to dC,,/dC, of —.111 was obtained at low 
angles and very good stability was indicated for the 
landing conditions. Completely removing the nacelles 
increased dC,,/dC;, to —.16, giving satisfactory sta- 
bility at all angles of attack. 

It was obvious that the engines and nacelles were 
the cause of the instability so an investigation of cowl- 
ings, nacelle location, shape and size was begun. Only 
slight improvement could be found. <A study of the 
flow between the nacelle and fuselage by means of 
smoke and threads attached to the model showed that 
the flow burbled at 8° angle of attack. This turbulence 
spoiled the flow over the tail almost up to the stall. 
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Wing Flaps - 45° 


Fic. 5. The Effect of the Engine 
Nacelles on Stability. ( Model 
Experiments were made using an auxiliary airfoil in 
this region to prevent burbling but after trying six 
different locations for the airfoil and obtaining only 
negligible improvement, tests with a fixed airfoil were 
abandoned. Considerable increases in drag were noted 
with all fixed auxiliary airfoils in this critical location. 
It was thought that if an auxiliary airfoil were so 
arranged that it could float about a spanwise axis at 
minimum drag until the airplane reached an angle ot 
attack of 8° (where it would encounter a stop) that 
the drag would be only slightly increased. By investi- 
gation of three airfoil positions and nine limiting angu- 
lar positions, the moment curve shown Test 15 
on Fig. 5 was obtained but with a drag increase of 
8% for the ship. 
Three type of windshields were 


as 


investigated but no 
changes noted in stability. 

Five different wing fillets were 
very large to small in size but for none of the cases 
tried was the stability as good as for the case with no 
fillet. The wing-fuselage junction on the Electra makes 
an angle of almost 90° so only a small fillet should be 
necessary. None of the fillets used increased the slope 
of the lift curve or the maximum lift a noticeable 
amount. From smoke tests, it seemed that the stream- 
lined flow did not have time to reestablish itself over 
the fillet after being disturbed by the nacelles. The 
fillets then, merely directed the turbulent flow on to 
the horizontal tail and it so happened in the case of 
the Electra that the tail was so located that with the 
fillets off, the wake missed the tail in the best manner. 
The above does not apply in general to a single engine 


tried ranging from 


airplane. 

When the trailing edge wing flaps were attached to 
the model, the slope of the moment curve was materially 
reduced at all angles of attack because of the increased 
downwash and lower tail efficiency. 

It occurred to the writer that the use of endplates 
on the horizontal tail would improve the stability by 
increasing the effective aspect ratio of the surface. As 
finally developed, these endplates were made to serve 
as fins and rudders so that the original vertical surface 
was removed entirely. The final moment curve 
shown in Fig. 8 compared with the original curve. 


is 


A 


Fic. 6. Effect of Fillets on Stability. 


Fic. 7. Effect of Wing Flaps on 


Test. ) Stability from Wind Tunnel Tests. 


substantial increase in the slope of the tail lift curve 
gave a good increase in stability at all angles of attack. 
The final value of dC,,/dC, in the cruising range was 
—.13 which was considered fairly adequate at that 
time. 

FiicHt Tests 

The problem of correlating wind tunnel and flight 
test data is one of considerable difficulty because of the 
number of variables involved. Inasmuch as the value 
of dC,,/dC;, had been used as stability criterion in 
the tunnel, it was decided to use the same criterion for 
flight testing. After several flights it was found that 
dC,,/dC, could be derived for any given C.G. posi- 
tion and throttle setting in the following manner. With 
some given C.G. position, say 34%, the ship would 
be trimmed for level flight. Pulling the plane up to 
some lower climbing speed, the elevator angle, stick 
force and airspeed would be read in steady flight. This 
was done for various climbing and diving speeds. The 
observer would then move the ballast in the cabin to 
some new position so that a 32% C.G. was obtained. 

With the original elevator tab setting, the test pro- 
cedure would be repeated for the same throttle setting 
and speeds. From the above data, the value of 
dC,,/dC;, could be derived. At an average C.G. posi- 
tion of 33% for any given velocity, AM (caused by 
moving ballast) = IV x d. 

dC, = AM/ Age. 

Plotting the elevator angles against velocity, at the 
given speed, d@ (change in elevator angle) is easily 
obtainable. The value of d@/dC;, can also be obtained 
by plotting the elevator angle against lift coefficient. 

Then dC,,/dé@ xX d0/dC, dC,,/dC;,, at the given 
speed, C.G. position and throttle setting. In_ this 
manner a direct comparison is obtained to the wind 
tunnel test. 

The wind tunnel value of dC,,/dC, of —.13 in the 
cruising range was found to be —.03 at the same value 
of Cr in flight. 


Elevator angles were measured by a pointer attached 
to the outer end of the elevator and noting the angle 
Using opera 


on a protractor scale painted on a fin. 
glasses, the angle could be read from the cabin to 
within .3 of a degree. Stick forces measured with a 
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i 34% CoGe 
Pinal 


Fic. 8. Results of the Wind Tunnel Study 
on Stability of the Lockheed Electra. 


spring scale are correct within 14 of a pound at high 
speed and '4 pound at low speeds. Power off read- 
ings were not taken at high speeds because the engines 
turned up too fast with closed throttles in the dive. 

Having taken stick force data, an oscillation test 
was run by pushing the stick ahead suddenly and releas- 
ing it. In this manner, an elevator-free oscillation was 
developed which was accepted as the final check on 
stability. High speed airplanes generally make very 
long oscillation even in their most stable loading condi- 
tions. Fig. 10 shows the stick force curve for the 
Electra with full load. If the stick force is divided by 
the square of the velocity, a curve is derived somewhat 
comparable to the tunnel test curve in that the forces 
are all expressed for a constant velocity. The Army 
uses this type of curve as an expression of the basic 
stability of an airplane. 

The curve showing the approximate flight path is only 
typical, as the amplitude of the oscillation is dependent 
on the amount of the disturbance which the pilot gives 
the control column, 

Little difference in time of oscillation 
approximately by use of a stop-watch) could be noted 
between conditions of locked and free elevators due per- 
haps to the fact that the elevators have such a high 
degree of static balance. A more unbalanced elevator 
would undoubtedly aid free elevator stability but the 
Department of Commerce requirement of a dynamic 
balance coefficient of .O8 determined largely the amount 
of counterbalance. 

As predicted by the wind tunnel tests, no trouble 
occurred in obtaining stability at high speed. Insta- 
bility with power on occurred with a 33% C.G. loca- 
tion at airspeeds around 90 m.p.h. with either free or 
locked controls. In this condition a push on the stick 
of one pound was enough to maintain a steady climb, 
but if the force was not applied, the ship stalled. A 
study of the airflow at the wing-fuselage junction by 
Inasmuch as 


(measured 


means of strings showed a good flow. 
the ship was stable power-off at this airspeed, the insta- 
bility must have been due to the action of the slip- 
stream in reducing the angle of attack of the tail as 
shown in Fig. 3. 


Fic. 9. Finding d0/dV_ Fic. 10. 
for 33% C.G. 


Full Load Stick Forces and an 
Approximate Elevator-Free Oscillation 
Path. 


It was thought that a large wing fillet might effect 
stability beneficially with the propeller acting although 
tunnel tests did not show them advisable. By actual 
test, the stability was slightly better with fillets off than 
on. Throughout most of the flight range, there was 
no change in the stick forces from the case with no 
fillets but at high angles, an increase of the unstable 
stick force of one-half pound was noted with the fillets 
on. The above should not be generalizéd, as a more 
thorough study of the problem may bring about op- 
posite results by the use of a different fillet. 

Changes made in the ship in the line of equipment 
and furnishings had brought the C.G. back to 34%. 
After a good deal of flight testing, it became obvious 
that the bump in the stability curve could not be elim- 
inated by simple changes. The easiest solution to the 
problem was found by increasing the tail span slightly. 
The allowable rearward C.G. position for stability 
ended up at 34.3% which is quite good for a multi- 
motored ship. 

The stability with flaps down and gear up is about 
neutral for the above C.G. with power on. With power 
off or with the wheels and flaps down, stability is good. 
The loaded C.G. position moves ahead over 1% when 
the landing gear is lowered on the Electra. (Fig. 12.) 

The longitudinal stability question cannot be settled 
by merely moving the C.G. ahead on the M.A.C. with- 
out running into trouble on landing with full forward 
loading conditions. Making a landing with a 16% C.G. 
with flaps fully down, the C.P. of the wing is so far 
behind the C.G. that the elevators cannot develop 
enough down load to force the tail down for a three 
point landing. With a low wing airplane especially, 
the downwash which aids in forcing the tail down 
is missing near the ground, so in case of a ship with 
a fixed stabilizer practically the whole burden of getting 
the tail down is thrown on to the elevators. In the 
case of the forward loading condition for the Electra, 
the ship could be stalled in the air (flaps and gear 
down) with an elevator angle of 13°. To make a three 
point landing in the same condition, required an ele- 
vator angle of 29°, the difference being due to the 
effect of the ground. With an elevator angle of 32° 
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38.34 CG. 


7 Power off 
21. 6% CoG. 


Indicated Speed mpeh. 
140 100 


34e3% CoG. 


indicated airspeed - mph. 
160 120” 
i 


i 


Variation of Stick Forces Fic. 12. 


with C.G. Position. 


Fic. 


Effect of Lowering the 
Landing Gear and Wing Flaps 


Fic. 13. 


with a Constant Loading | 
Condition. 


the tail stalled. A careful investigation of the best 
limiting elevator angle for landing is well worth while. 
It is important to make a three point landing in all 
loading conditions to minimize the chance of nosing 
over on a muddy field or on sudden application of the 
brakes. 

The arrangement of the tail surfaces on the Electra 
is fortunate in that a section of the upper surface of 
the tail cone can be hinged as shown in Fig. 14 to act as 
up-elevator. No spring device is used to keep this 
section closed in flight so there is no sudden increase 
in stick force when the elevators pick up this center 
section on the up throw. By closing the gap between 
the elevators in this manner, the longitudinal control 
becomes very effective on landing so that a greater 
range of C.G. travel is allowed. 


CONCLUSION 
The author regrets not having had time to work up 
all the data in connection with the effects of flaps and 
throttle setting on dC,,/dC, but before doing so he 
would like to develop more accurate means of obtain- 
ing data on elevator angles, flight path in oscillations 
and so on. He would like to point out one factor 


which has generally been overlooked in connection with 
stability. An airplane is not satisfactory if it is mar- 
ginally stable or slightly unstable in a climb at the 
low altitudes, for the indicated airspeed in a clim) at 
this altitude becomes the cruising indicated airspeeds 
for high altitude flying with a modern airliner. The 
pilot is then forced to fly a plane which hunts in cruis- 
ing as the longitudinal stability is a function of the 
indicated airspeed. 


Fic. 14. Tail Cone Flap on the Electra. 
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TUDENTS of fluid mechanics agree now rather 

generally that the separation of the boundary layer 
from a wing surface is the cause of burbling and is re- 
sponsible, therefore, for the complicated phenomena near 
the stall. There is consequently a great practical as 
well as scientific interest in obtaining an analytical 
method of predicting such separation. One method 
of analyzing flow separation in the case of a turbulent 
boundary layer was proposed by Gruschwitz' and 
seemed to offer a practical solution. The limited 
experimental evidence in favor of this method, how- 
ever, was not sufficient to command confidence in its 
general application. In the present research an 
attempt was made to verify the Gruschwitz method by 
an extended experimental examination of the turbulent 
boundary layer over a typical wing profile at high 
Reynolds Numbers. The data here presented lead 
to the conclusion that the Gruschwitz method does not 
describe the actual phenomena and is, therefore, not 
the general solution hoped for. The evidence now 
seems to point to a dependence of flow separation 
upon several other variables besides the pressure 
gradient and the thickness of the boundary layer which 
need further experimental examination. Such investi- 
gation is now under way at the Massachusetts Insti- 
tute of Technology and a further report will be made 
in a subsequent paper. 

The mechanics of flow separation depends upon the 
motion of the air within the boundary layer. The 
particles within the boundary layer are retarded by 
the surface friction and are usually carried agains: a 
rising pressure by the momentum transferred to them 
from outside of the boundary. If the pressure rise 
is sufficient the particles lose energy, come to rest, and 
separation occurs. If the momentum is transferred 
from the exterior stream by molecular exchange 
purely, the motion is defined as laminar, while the 
motion is defined as turbulent if the exchange is molar. 

All boundary layers begin at a stagnation point and 
remain laminar for a certain distance downstream 
whatever may be the character of the exterior flow or 
of the surface. If the prevailing velocity is sufficiently 
high and the surface sufficiently large, however, the 
layer will become turbulent and continue so throughout 
the remainder of its length. The chief characteristic 
which distinguishes the turbulent flow from the laminar 
is its much greater ability to transport momentum 
laterally. This gives it a correspondingly greater 
ability to proceed against a rising pressure. 


1E. Gruschwitz, Ingenieur Archiv., Vol. 2, 1931. 
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The separation of laminar boundary layers has 
recently been studied by von Karman and Millikan? 
who explained the variation of the maximum lift 
coefficient with low Reynolds numbers in the case of 
certain wing profiles. The separation of the turbulent 
boundary layer, however, is of greater importance at 
the higher Reynolds numbers used in engineering since 
it is responsible for most of the burbling phenomena 
met in practice. 

The studies presented here are an outgrowth of a 
recent investigation by Gruschwitz* on the behavior of 
turbulent boundary layers for various pressure gradi- 
ents. Since this method is used for comparison with 
the experimental data here reported, it will be well 
to review it briefly. 


GRUSCHWITz’s APPROXIMATION FOR THE CALCULATION 
OF A TURBULENT BouNDARY LAYER 
The boundary layer has to satisfy the momentum 
law, first applied to the boundary layer by von Karman.‘ 
This equation is 


1 ad ld 6 1 dp To 


where 
4= distance from the stagnation point measured 
along the surface. 


4’ = distance normal to the surface. 
z= the velocity at any point in the boundary layer. 
V’ = the velocity outside of the boundary layer. 


6 = thickness of the boundary layer, for 
= %,(V —v)> 0. 


G, 


p= static pressure along the surface, (varies only 
with x and equals the pressure of the potential flow). 


p = mass density, t = shearing stress at y = 0. 
Eq. (1) can be rewritten with 
dz ds’ 


= 3, and (V —1) dy = 5° 


and gives then 


1d 1 @ To_ 
Introducing 
= (l dy 


gives the final equation 


2C. B. Millikan, 4.S.M.E. Transactions, 1934. 
%E. Gruschwitz, /ngenicur Archiv., Vol 2, 1931. 
4Th. von Karman, Z.4.M.M., Vol. 1, 1921. 
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2 
@ has the dimension of a length and can be used as a 
measure for the thickness of the boundary layer. It 
is proportional to the loss of momentum since pV’*6 is 
the difference between the momentum of the mass 
inside of the boundary layer flowing with the potential 
velocity (lV), pV and the actual momentum 

Eq. (2) is essentially a differential equation for 4 
and can be solved for a known potential flow along the 
surface if reasonable values are assumed for 1./pV* 
and 3*/8. 

To obtain information about the velocity distribu- 
tion in the boundary layer Gruschwitz assumes that 
the curves v/V = f (y/#) at various points, +, belong 
to a family of curves with only one parameter. In 
this case v/V at any fixed level y/@ might be used 
as the parameter. It proved to be convenient, however, 
to choose the parameter 4 = 1 — (zv:/V)°* evaluated at 
the level y = 4 instead of the velocity ratio v./V’ itself. 
The experimental values of 3*/@ as determined from tur- 
bulent boundary layer test plot as a single curve against 
the parameter 4. This relates the center of gravity to 
the area /* (l1—v/V)dy with the parameter of the 
velocity distribution curve, since 


V) (1 + v/V)dy 
1) 6 
j (1—v/V) dy 


This relation between the center of gravity and the 
parameter can be taken as a substantiation of the 
assumption of a one parameter family of curves. 

Considering now a particle with the velocity, w, at 
a distance y = 4 from the surface Gruschwitz concludes 
that the variation of its energy with respect to x is a 
function of 7, V, 9, p,and of the viscosity v 


d/dz (> = = f (n,V,0,0,») 


By dimension analysis the following relation can be 
6 dg 


derived 
Ve Ve 
V2p/2 dz ) (> ) @) 


When measured values of the terms in Eq. (3) were 
plotted, the influence of ’ 6 /y was found to be negligible 
in the range tested and (26/V"p) (dg:/dx) to be a 
linear function of 4, that is equal to a 4 + Bb for all 
values of .r upstream of the region of separation. Eq. 
(3) can then be written 


V2p/2 dr (4) 

where a and b were found to have the numerical values 
a = 0.00894 and b = — 0.00461. 


Having the empirical differential Eq. (4) for 4 which 
determines the form of the velocity distribution curve 
and the Eq. (2) for 4 which determines the boundary 


Hi. PETERS 


layer thickness, the two equations can be solved simul- 
taneously to yield the approximation of the turbulent 
boundary layer if the pressure distribution along the 
surface is known. The solution is reached by succes- 
sive approximations. First, Eq. (2) is solved for 6, 
then with 6 known Eq. (4) is solved for 7. To make 
the first solution of Eq. (2) possible, assumed values 
of @*/@ and 7,/pl~ have to be used. One notes that 
To/p V* is a relatively small quantity compared with the 
other two terms of Eq. (2) provided the order of 
magnitude of dV*/dx is appreciable. The assumption 
that 7,/pl’* is constant, the order of magnitude of which 
is estimated, will therefore not introduce a large error, 
2*/6 also can be assumed to be constant for the first 
approximations. The second approximation will be 
closer since 8*/@ is approximately known as a func- 
tion of 7. In addition to better values for 8*/@, the 
estimated values 7,/pV* can be corrected by using the 
known values for a flat plate. 

The integration constant, or the 4 value at the 
beginning of the turbulent boundary layer, is approxi- 
mated by the assumption that the thickness of the 
turbulent boundary layer is equal to the thickness of 
the laminar boundary layer at the transition point. The 
laminar boundary layer might be calculated with the 
Pohlhausen® method. The value for 4 at the transi- 
tion point is of minor importance since the 4 curve 
converges very rapidly. Gruschwitz suggests the value 
4 == 0.1 at this point as a good assumption. 

The separation of the turbulent boundary layer was 
assumed to occur when the velocity ratio v/V at a 
fixed level in the boundary layer reached a minimum 
value. Gruschwitz concluded from experimental data 
that the separation with turbulent boundary layers 
occurs with v:/V = 0.447 at the level y = 4 correspond- 
ing to O.8. 


M. I. T. Tests anp RESULTS 

The experiments to test the validity of the Gru- 
schwitz method, especially at higher Reynolds numbers, 
were made at the Massachusetts Institute of Technology 
on a symmetrical wing section of 7% ft. chord whose 
span extended from wall to wall of the 7% ft. wind 
tunnel. The maximum thickness of the section was 
1/6 chord lengths and occurred at 33% of the chord. 
The general arrangement can be seen in Fig. 1. The 
surface of the side of the wing on which the measure- 
ments were made consisted of thin brass sheets, care- 
fully fitted together to give a continuous smooth sur- 
face. The large size wing section changed the normal 
velocity distribution in the wind tunnel completely 
and therefore made the determination of Vx, the velo- 
city of approach, quite doubtful. 

To study the influence of various pressure gradients 
the boundary layer was measured on one side of the 


5K. Pohlhausen, Z.4.17.M., Vol. 1, 1921. 
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Fic. 1. General Arrangement of the Wing Section in the Wind Tunnel. 


to +11 
The measurements were taken in the middle section 
flow could be 


wing at angles of attack varying from —9 


of the wing where two-dimensional 
assumed to exist. 

The static pressure distribution was measured by 
means of 36 holes of 0.8 mm. diameter located at 
intervals along the surface. Fig. 2 shows p/g, = f 
(*%) giving the pressure distribution for various 
angles of attack plotted against the surface distance 
r% = (x/s) 100 where x is the distance from the 
leading edge, S the total distance from leading edge to 
the trailing edge, and gp = p+ the total pres- 
sure outside the boundary layer. Provision was made 
for the mounting of a pitot tube at each of 13 stations 
along the centre line of the wing and for its accurate 
setting at any desired level from the surface with the 
precision of + 0.02 mm. by means of a micrometer 


Fic. 2. Pressure Distribution for Various Angles of 
Attack. 


screw. The smallest distance possible, namely, the 
distance equal to the radius of the pitot tube, was 
determined by the closing of an electrical circuit at 
the moment that the tip of the pitot tube touched the 
surface. The pressure readings were assumed to equal 
the pressure prevailing at the centre of the opening 
of the tube. Various sizes of tubes were used; their 
readings compared quite satisfactorily with each other. 
The smallest pitot tube employed had an outside dia- 
meter of 0.46 mm. 

The procedure of the test was to measure first the 
pressure distribution along the surface for a fixed 
angle of attack for two or three wind velocities and 
then to measure the total pressure, g, at various levels 
from the surface at each of the 13 stations, one sta- 
tion at a time. Repeated checks on both the surface 
pressure and the total pressure distribution across the 
boundary laver proved the consistency of the flow. 

From the two equations whose right sides are known 
from the measurements, vp/2 = g—p for the flow in- 
side the boundary layer and l’*p/2 = g,—p for the poten- 
tial flow just outside, the terms (v//”)* and v/lV are 
calculated and plotted against the normal distance y 
planimetered areas give 


from the surface. The 


— fi 1—v/V)dy anda = 1— (v/V)*)dy 


l'rom these two values the value of @ is obtained by sim- 
ple subtraction since A — f v/V (1—v/V) dy. 


The Gruschwitz parameter for the velocity curve, 
4—1—(vV)? is found by reading the value 
at y = 6 from the plots of (v//")* against y. The quan- 
title 8*/@ plotted against 7, can be represented by the 
same single curve that Gruschwitz obtained. 

The shearing stress at y = O in terms of pl’, can 
be obtained either by von Karman’s momentum Eq 


(2) or from 7,/pl” = (v/V*) (dv/dy), (y=90). Both 
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Fic. 4. Velocity Distribution in the Boundary Layer at Points Upstream of the Trailing Edge. 
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methods are applied and the results plotted against the 
Reynolds numbers Ro = V6/v in Fig. (3). The com- 
parison of these values of the shearing stress with 
those for flat plate shows clearly that the shearing 
stress is less for surfaces with a positive dp/dx than 
with dp/dx = 0, although the scattering of the results 
is fairly large due to the inaccuracy of the determination 
of dé/dx in the one case and of dv/dy at y=O in 
the other. 

The transition from laminar to turbulent boundary 
layer takes place quite consistently at = 600 to 650; 
the transition can be clearly identified also from the 
change of the 4 values against +% (see Figs. 5 and 6) 
as well as from the form of the velocity distribution 


curves themselves. 
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Fic. 5. Comparison between Measured and Calculated 
Values for the Boundary Layer Thickness 9 and 
the Velocity Parameter 7. 


Fic. 6. Comparison between Measured and Calculated 
Values for the Boundary Layer Thickness @ and 
the Velocity Parameter 1. 


The point at which the turbulent flow separates from 
the surface is determined in the manner indicated in 
Fig. 4 which shows the velocity distribution near the 
surface at a few stations upstream of the trailing edge 
when the angle of attack had the values of 2 = 5°, 7°, 
9° and 11°. There is some uncertainty about the velo- 
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city values measured near the separation point due to 
the unsteadiness of the flow in this region. The separa- 
tion may occur at a surface position somewhat further 
upstream but certainly not downstream of the point 
determined by Fig. 4. 

To compare the values of @ and 7 as obtained from 
the measurements with the values computed by the 
Gruschwitz method these two values were calculated 
from the two differential Eqs. (2) and (4) for the 
surface pressure distribution as given in Fig. 2. The 
full-drawn lines in Figs. 5 and 6 represent the values 
thus obtaintd for 2 = — 9°, and + 9° by the first and 
second approximation. The plots for «=— 5°, 0°, 
+ 5° and +7° are omitted here, since they show 
essentially the same agreement for the measured and 
computed values as the plots for a =-—9°. The 
marked points indicate the experimental values evalu- 
ated from the measured data. To solve the two differ- 
ential equations the Czuber’s® graphical method was 
used. This method is based on the fact that all line 
elements for += const. of a differential equation 


dy/dx + yfi(x) = fe(x) 
lead to one point whose coordinates are 
t+ l/f,(@) and Y=f,(2)/f,(2). 


For the first approximation the shearing stress and 
the ratio 4*/@ were assumed to be constant for the 
whole surface having the numerical values 

t>/pV* = 0.0017 and #*/6=—1.4. For the second 
approximation the values=f (7%) and the flat plate 
values 7,/pl* = f (R@) were used. At the transition 
point 8 was assumed to equal its measured value at 
that point and the 7 curve has been started with 4 = Y. 

The comparison between computed and tested values 
of 4 and @ showed that for the range of angles of 
attack, for which no flow separation occurred, a fairly 
good agreement was found (see Fig 5 for example). 
Fig. 6, for the angle of attack of + 9°, on the other 
hand, shows that when separation occurs the calculated 
values of 4 and @ depart more and more from the 
measured values as the region of separation is 
approached. Separation is seen to occur approximately 
at x % = 75, and with 7 = 0.8. It is noteworthy that 
the value of 4 as calculated does not reach the value 
4,=0.8. Therefore, one cannot predict from the cal- 
culation anything in regard to flow separation since, 
according to Gruschwitz, only 4 = 0.8 would indicate 
the danger of separation. Similar results were obtained 
fora=+11°. 

The essential reason for the disagreement between 
the computed and the measured values of 7 and @ is 
found in the validity of the linear interpolation which 
is made in setting up the differential Eq. (4) for the 
form parameter, 4, of the velocity distribution. Fig 7, 


°E. Czuber, Z.S.f. Math. and Physik, Vol. 44, 1899. 
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where (26/l’"p)dg,/dx is plotted against 7, shows that 1 
the functional relation is not a linear one, a 4+), SI. 

a =-9° 
even for cases of no separation except in a first a: ace 
rough approximation, and that this approximation 
is not valid near the region of separation and > ae 
for some distance upstream. Furthermore, the [ + a= 

| 


range of 4 for which the linear interpolation is roughly 
valid is uncertain due to the fact that the separation 001 | 
point seems to jump upstream rather than to move 
gradually with increasing angle of attack. 

It appears therefore, from these studies that the 
Gruschwitz method cannot be-used to determine the 


location of the separation point or even in many cases 


to predict whether separation will occur at all. There 
is need for further experimental information concern- 
ing the relation between the velocity distribution in 
the boundary layer, pressure gradient, thickness of 


04 5 
the boundary layer, curvature of the surface and ; 
; Fic. 7. Experimental Data Compared with the Straight 

degree of turbulence of the approach flow. Line Interpolation as Suggested by Gruschwitz. 
4 
3 


| 
| | 
| | ‘ 
aN 
| 
| 0.00894 n~-0.00461 
| e| | 
| 
| ex | 
| 
\ 
| \ 
| | \ | 
| | | 
| 
| | 
(BY 
I 
2 
REEL 


rht 


SEPTEMBER, 1935 


JOURNAL OF THE AERONAUTICAL SCIENCES 


The Strength of Transverse Frames 
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(Received January 11, 1935) 


SUMMARY 

HE main transverse frames of many aircraft struc- 

tures are extremely redundant. In the hulls of rigid 
airships and flying boats and in the fuselages of air- 
planes the tendency of design is towards deep, con- 
tinuous ring frames, which are in themselves struc- 
turally redundant, and may present quite difficult prob- 
lems in stress analysis. The problem becomes much 
more complicated when transverse bracing members 
are added to the frame. 

A method of stress analysis of such structures is 
presented in this paper. The analysis of a main trans- 
verse frame of a rigid airship, constructed as a deep 
ring with tight radial wires, is given as an illustrative 
example. 


Continuous RiInGc FRAME WITHOUT CRross-BRACING 


A ring of any shape and with any loading can be 
made statically determinate by cutting it through at 
some point. <A stress analysis by the principle of least 
work or other strain energy method can be made by 
applying a moment, an axial force, and a shear at the 
cut section. If the structure and loading are symmetri- 
cal about an axis, and if the ring is cut at one of its 
intersections with that axis, there is no shear at the 
cut section, and the stress analysis requires the solu- 
tion of only two unknowns, the bending moment and 
the axial force at the cut. 


ConTINUOUS FRAME WITH Cross-BRACING 


Continuous frames with cross-bracing are usually 
treated in stress analyses as though they were pin- 
jointed at the intersections of the cross-bracing with 
the ring. If the ring is itself quite flexible, and 
dependent mainly upon the cross-bracing for its 
strength and stability, the assumption of pin-joints will 
probably give a sufficiently close approximation to the 
stress distribution. 

When the ring is deep and rigid enough to provide 
a considerable part of the total rigidity of the assembly 
of ring and cross-bracing, the pin-jointed assumption 
is no longer satisfactory.2, A frame of this type may 
be made statically determinate by omitting the cross- 


_1C. P. Burgess, Airship Design, page 203, The Ronald Press 
Company, 1929. 

® Miss H. M. Lyon, The Strength of Transverse Frames of 
Rigid Airships, Journal of the Royal Aeronautical Society, 
page 497, June 1930. 


bracing, and cutting the ring at a section, preferably 
on the axis of symmetry if there is one.* 


Fic. 1. Deep Ring Frame with Tight Wires. 


The main frame of a rigid airship, shown in the 
accompanying figure, may be taken as an example of 
a continuous deep ring frame with cross-bracing. It 
consists of a uniform, trussed ring, reinforced with 
tight radial wires. 

The loading system consists of a vertically down- 
ward load of 2 lb. (1 Ib. for each side) at the bottom 
of the frame, and tangential lift forces applied at all 
the outer joints of the frame, perpendicularly to the 
radius through the joint, and proportional in intensity 
to sin 8, where @ is the angle between the radius through 
the joint and the upper vertical radius. In a circular 
frame with » uniformly spaced joints, and carrying a 
dead-load of G Ib., the external tangential force at any 
joint is equal to (2G sin @) /n, or .111 sin @ Ib. in the 
case of the 36-sided frame with a 2 lb. load. 

The stress distribution is calculated by the method 
of least work. Two of the unknowns are the bending 
moment and the axial force at the top. The remain- 
ing unknowns are the stresses in each pair of the 
omitted radial wires. 


3E. E. Lundquist and W. F. Burke, General Equations for 
the Stress Analysis of Rings, N.A.C.A. Technical Report 509, 
1934. 
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HARACTERISTI F THE FRAM 

OVERALL RADIUS 106 FT. 
RADIUS TO NEUTRAL AXIS 100 FT. 
INSIDE RADIUS 88 FT. 
MOMENT OF INERTIA OF CROSS-SECTION 

OF RING 232 FT? iN? 
CROSS-SECTIONAL AREA OF WIRE .0235 In® 
E, FRAME 1.0 
E, WIRES 2.8 


zo 27 ©.03| 11.36 

30 13.40] 34.20 

40 1.50 23.40} 50.00) 11.36 

50 2.82 35.12| 64.28) 34.20 

4.87) * 50.00} 16.60} 50.00) 17.56 

T° 1.867) ©5.80| 86.60) 64.28) 34.20 

80 12.00) * 82.64) 93.97) 76.60) 50.00) 17.36) 

90 11.46; 100.00} 9848) 86,60) 64.28| 34.20 

100 | 2448) 111,36] 100,00} 93.97) 7660) 50.00) 

120 5 

130 | 55.54) 164.28| 86.60/ 98.48 86.60 

140 | 69.18; 11@.60| 76.60} 93,97} 100.00) 93.97) 16.60] 50.00) 17.36 

150 | 84.24) 186.60] 64.28) 66.60) 2848) 298.48) 86.60) 6478) 34.20 

160 | 100.48; « {93.91| 50,00} 76.60} 93.97| 100.00] 93.97) 76.60) 50.00) [7.56 
(70 | 117.55) 198.48) 34.20) 64.28) 86.60) 98.48|) 9846) 86.60) 64728/ 34,20 
160 | 13496; 200.00| 34.20] 64.28] 86.60) 98.48) 96.48) 8660) 6428) 34.20 


It is assumed that the wires have sufficient initial 
tension to take compressive forces without going slack. 


NoTATION 


In the frame shown in the figure, let 


R = radius of the ring along its neutral axis. 

I =moment of inertia of the cross-section of the 
ring. 

A =cross-sectional area of a radial wire. 

6 angle between the radius to any point and the 
upper vertical radius. 

a == value of 6 where a radial wire is attached to the 
ring. 

M. = bending moment at any section of the statically 
determinate (cut) ring. 

AM, = bending moment due a moment numerically equal 


to R applied at the cut section, 

M. = bending moment due to unit axial tension at 
the cut section. 

M, = bending moment due to unit tension in a redund- 
ant wire. 


M = bending moment at any section of the complete 
frame. 

Xx, bending moment at the top of the complete frame. 

X. == tension at top of the frame. 


Xs = stress in wire I. 


m,M, L/EI 4833 
L/EL 8332 
M,M,L/EL 3071 
MM, L/EI 3913 
MoM, L/EI 4288 
MoM, L/EI 
3574 
L/EI 2649 
MoM, L/EL 1578 
G07 
me 13,527 
mM, M, L/EL 13,499 
M,M, L/EL 8,586 
M, 8,068 
M,MsL/EI 1,093 
M,M,L/ET 5,781 
M, M, L/EL 4,269 
M,M, 2,798 
M, M, |,488 
MM, L/EL 


mM? L/eL 20,290 
9,760 
L/EI 10,940 
M,M,L/EI 10,872 
L/EI 9,681 
M.M,L/EL 1,665 
M,M,L/EL 5,242 
M,M, L/EI 2,88) 


M,M,L/EL 


S,'L/ea 
MyM, L/EL 6,395 S35,L/E€A 1,099 
M3M5L/EI 5,411 33 Bie 
M3M, L/EI 4,104 $35,L/EA 
M3sM,L/EI 2,147 $35,/EA ° 
MsMoL/EI 101 Ss - B16 
MsMoL/EI Sy -1,099 
6,619 S,* 2,217 
6,069 L/EA 117 
MaM, L/EL 4,124 S454 L/EA 362 
3,586 L/EA ° 
M4M,/EL 2,223 382 
MaM5L/El 1,099 S4 Se L/EA 
M4M,L/EL 348 - 966 
L/EL 5,995 St L/EA 1,783 
MsMcL /EI 5,221 L/EA 283 
3,993 $554 ° 
2,013 - 283 
MsMgL/EL 1,364 L/EA - 533 
mM 4,85! 1,401 
MoM, L/EL 3,918 $3.5, /EA ° 
MgMg L/EI 2,688 S_L/EA - 16! 
1,465 S_ Soi - 283 
3, 5,L/EA - 362 
M4? L/EL 3,370 L/ea. 1,250 
L/EL 2,439 ° 
/EL $02 $,5,,L/EA ° 
M,? 1,895 5," L/EA 1,401 
L/EL 1,146 282 
M,M,U/El 454 SU/EA 382 
| Mgt 164 /ea 1,783 
| /EX 33 =| SsS_L/EA 118 
| iss] 2,217 
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X. = stress in 
Xs =stress in 
X. = stress in 
X; = stress in 
Xs = stress in 
X» = stress in 
Xw = stress in 


wire II. 
wire ITI. 
wire IV. 
wire V. 
wire VI. 


THE 


wire VII. 
wire VIII. 


STRENGTH 


OF 


The characteristics of the frame are given in Table 1. 


M, R. 


M: = R (1—cos @). 
M,=R (sin@—sina), 
Where @ is less than 2, M, = 0. 


CALCULATION OF THE STRAIN ENERGY TERMS 


The strain energy terms for the least work equations 


are of the form 


great. 


f (M,, M,,/El) dex 


The labor of performing all the integrations for the 
eleven bending moments, 1/, to M,, would be very 


Good approximations to the integrations may 


be obtained quite simply by dividing the ring into 
arcs of 10°, in each of which a constant bending 
moment, equal to the mean for the arc, is assumed. 


The strain energy term for each arc may then be 


expressed in the form 
M, M,L/EI 


TRANSVERSE FRAMES 


where L is the length of the arc along the neutral axis. 
The bending moments are tabulated in Table 2. 
The strain energy terms of the form =M,,M,L/EI are 
calculated by multiplying the bending moments in the 
horizontal rows of Table 2 together in pairs, adding the 
resulting columns, and finally multiplying the additions 
The values of 3M,,M,L/EI 
computed in this manner are summarized in Table 3. 
The strain energy in the wires must be computed in 


by the constant L/E/. 


addition to the work of bending the ring. 


The terms 


for the wires are of the form S, S, L/EA, computed 
from the stresses, S, given in Table 4. The summa- 
tions of the energy terms are included in Table 3. 


Ss 


Se 


NO. x, Xe X3 X4 Xs X, Xs Xe CONSTANT CHECK 
\ 13,527 13,499 8,587 | 8068 | 7,093 | 578! 4,289 | 2,198 1,488 SI@ 4,833 10,419 
é '3, 499 20,290 9,160 | 10,940 | 10872 9,68! 1,065 | 5,242 2,881 1,020 8, 332 100, 182 
3 8,587 9,760 9,285 | 1,494 | 6,227 | 4,538 | 2,147 1,13) J- 11S 902 3,071 S\, 823 
8,068 10,940 1,494 | 8836] 6,786 | 5,106 3,586 1,841 382 o18 3,913 56,334 
5 1,093 10,872 6,227 | 6,786 | 17,778 | 5,504 | 3,993 | 2,330 83) 262 4,288 55,440 
6 5,781 9,681 4,538 | 5,106 | 5504] 6252 | 3918 | 2,537 1,182 128 4,\e1 48,188 
1 4,289 1,665 2,147 | 3,586 | 3,993 | 3,918 | 4,620 | 2,439 1,389 $02 3,514 38,722 
8 2,798 5,242 1,131 1,841 2,330 | 2,537 | 2,439 | 3,296 1,429 Ble 2,649 26, 508 
3 1488 2,881 115 382 83! 1,182 1,389 | 1.429 | 2,547 1,031 1,578 14,623 
516 020 |- 902 |- |- 262 128 502 BIG 1,03) 2,312 601 5,210 


10,333 
3a} 63480 8,569 5,45) §,122 4,503 3,670 2,123 1,176 945 328 3,068 44,140 
42} 59644 6,051 5,122 4,812 | 4,231 3,448 2,558 1,669 888 308 2,883 42,036 
Sa} 52436 1,018 4,503 4,231 3,119 3,031 2,249 1,467 180 27 2,534 36, 956 
Ga} 42137 5,169 3,670 3,448 3,031 2,411 1,833 1,196 636 221 2,06S 30, 121 
Ta} .31107 4,280 2,123 | 2,558 2,249 1,833 1,360 887 412 164 1,532 22, 341 
8a} .20685 2,192 i, iT | 1,669 1,461 1,196 887 519 308 107 1,000 14, 519 
9a! 11000 1,485 345 180 C36 412 308 532 7, 153 


2,689 


29,850 

3834 | 2312 | 1,124 868 24 645 |- 1,060 |- 1,230 3 1,081 
13 2,889 2,312 | 4,024 | 2555 | 1,658 | 1,028 112 + S506 926 1,030 14,296 
14 3,194 L124 | 255s | 4059 | 2473 | 1.1744 863 oT 533 L154 18,484 
3,912 ses _ | 1658 | 2473 78! | 2085 | 1.341 _546 23 2,096 18,667 
3,385 24 |, 144 2,085 | 3,260 1,552 217 338 2,042 16,3715 
2,450 |- 645 172 | 134: | 1,552 | 2717 1,421 7109 1,649 11,929 
1,396 |- toco Soe 51 546 217 1,121 2,383 274 1,046 6,868 
2,519 


5,214 


Sa A2367 505 | 1,224] 1,607 | 1,657 | 1,434 | 1,038 59! 214 1,487 12,641 
a 5S639 | 1,607 | 2,111 2,177 | 1,883 | 1,363 1717 281 1,952 16, 608 
Sa -57369 | 1057 | 2,177 | 2244 [ 1942 | 1.406 601 290 2,013 17,125 
lea 49641 591 1434 | 1,883 | 1,942 | 180 | 1,216 693 251 1,742 14, 818 
Te 35929 428 1038 | 1363 | 1,406 1.216 502 181 1, 261 10,725 
Ba -20412 244 591 1717 801 693 502 286 103 118 PA 
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The terms of the standard least work equations are 
the sums of the strain energy terms from the ring 
and wires, neglecting the work of the axial force and 
shear in the ring as being unimportant, although it is 
quite practicable to include them if a high degree of 
accuracy is required. 


The numerical least work equations for the ten 
unknowns are written in tabular form in Table 5. They 
may be solved with the aid of a calculating machine in 
a few hours. The method recommended is a modifica- 
tion of that invented by the astronomer Gauss.* As 
it appears to be very little known, a brief description 
of the procedure may be worth while. 

Starting with Eqs. (1) to (10) in Table 5, the first 
term in Eq. (1) is divided successively into the first 
terms of Eqs. (2) to (10), giving the fractions listed 
in the column preceding Eqs. (2a) to (10a) in the table. 
These fractions are then multiplied into all the terms 
but the first one of Eq. (1), and the products are the 
terms of Eqs. (2a) to (10a). Egs. (11) to (19) are 
obtained by subtracting the terms of Eqs. (2a) to (10a) 
from the corresponding terms of Eqs. (2) to (10). 


The procedure for reducing Eqs. (11) to (19) to 
Eqs. (20) to (27) is the same as in obtaining (11) 
to (19) from (1) to (10), and so on until there is 
but one unknown equated to a constant. 


An important part of the procedure is the check 
term which guards against an error being carried 
through the calculations. The check term in Eqs. (1) 
to (10) is the algebraic sum of all the other terms 
in the equation. The check term in Eqs. (11) to 
(19) is computed in the same way as the other terms, 
and if the work is correct, it will still be the sum of 
all the other terms in the equation. 


4C. P. Burgess, The Application of the Principle of Least 
Work to the Primary Stress Calculations of Space Frameworks, 
Transactions of the A.S.M.E., Aeronautical Engineer, Vol. 1, 
No. 3, page 107, July-September, 1929. 
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TABLE 6. 


BENDING MOMENTS AND AXIAL LOADS DUE TO 2 LBS. 


WITH WIRES WITHOUT WIRES 
m/R Pp u/R Pp 
LBs. LBs. LBS. Les. 

0 .0563 | -.4158 | .1563 | -.5146 
10 .0506 1491 
20 .0502 | -.3692 .1280 | -.4456 
30 .0398 .0944 
40 .0294 | -.2312 .0509 | -.2514 
50 .0129 0007 
60 |-.0070 |-.0154 | -.0523 
10 | -.0267 - 
80 | -.0525 .2505 | -.1490 
90 | -.0706 -.1835 
100 | -.0961 .5284 | -.2028 .6364 
110 | -.2029 
120 | -.1180 -7617 | -.8070 | .8345 
| -.1018 -.1337 
140 | -.0728 .9053 | -.0607 -8942 
150 | -.0178 .0385 
160 | .0726 | .1629 
170 . 1886 
180 . 3564 .6371 ATCT 5146 


The. solutions of the equations are: 


X, = 0563 X, = — .0569 
A, = AS X, =— .1104 
A, = 1D X, =— .1107 
X, = 0483 X, = — .0554 
X, =— Xy= 0475 


The tension in the bottom wire to give equilibrium 
with the other wires is .1488 Ib. 

The bending moments and axial forces of the deep 
ring frame, with and without the radia! wires are com- 
pared in Table 6. The wires constitute only about 
3% of the total weight of the frame; but they relieve 
it of about one-half of the mean bending moment in the 
unwired frame. With heavier wires or a frame hav- 
ing a smaller moment of inertia, the fraction of the 
bending moment taken by the wires would be greater. 
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The Boundary Layer and the Drag of a Body of Revolution at Large 
Reynold’s Number 


C. C. Feprarvsxy,* Central Aero-Hydrodynamical Institute, Moscow 
(Received November 1, 1934) 


HE Karman integral relation deduced for a two- 

dimensional flow will assume for the boundary 
layer of a body of revolution at a steady motion the 
following form: 


d 6, d 
putdy — 27 r’ pud 
AS, r’ puedy wil r’ pudy 


diJo 


where / is the length of the generator of the body of 
revolution measured from the bow point of the body. 
7” = radius of an element of the surface normal to 
the body of revolution (see Fig. 1). 
p = fluid density. 
y= distance from the surface of the body, meas- 


ured along the normal axis. \ 
= boundary layer thickness. 
u = velocity inside the boundary layer. 
us = velocity just outside the boundary layer. 
p = static pressure. 
r = radius of a body of revolution. — 
t= frictional force per unit of the surface area. rm % 
=r 8 = 8 =— w 9); 
u, = us (1); p =p (1) andr =r (1) We shall then have: 
let us introduce the following notation : df ldr 
V = velocity far from the body. af, 
=. static pressure corresponding to the velocity 
P—Po 
“F *~ tore In the case of a turbulent boundary layer let us 


denote the velocity distribution and the skin friction 


x= axial distance from the bow point. be the follow! : 
the by the following expressions. 
¥= coefficients of kinematic viscosity. “— w(4)’ (3) 
Re= Reynold’s number checked lengthwise on the : 
body. =K K pita (4) 
v v 


From Fig. 1 we can see that: 
where “nx” and “A” are functions of the Reynold’s 


Number and m = — 2n/1 + n. 
where 6 is the angle between the tangent to the gene- . 
U Eqs. (3 1 (4) Eq. (2 d 
rator of the body of revolution and the longitudinal 


to the following form: 


Assuming that even in the case of a maximum value + = | 28 (1—e) 
of the quantity “y” in the boundary layer ( i. e. &) it is 
small in comparison with “r’”’, Eq. (1) can be consider- =-—Kki-m(L)" FF 
y simplified by substituting “r’” for “r. 
where 
*C. C. Fediaevsky is an aerodynamical engineer engaged B= 1 1 


since 1925 in research work at the Central Aero-Hydrodynamical ~ ae 1 
Institute, Moscow, and also serving as Professor at the Airship = sl 
Institute. His paper was sent to the Journal by Prof. A. and 

Nekrassov, Vice Director of the C.A.H.I. “as a contribution 

to the problem mentioned by von Karman in his article in Om ae, 
the first issue of the Journal of the Aeronautical Sciences.” 2n+1 n+l 
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Integration will give the following results: 
1 m m l—e 
6 = R L F 
| ( ‘) 


(—e) (1—m) 


m 
2 2B 
x| fir 
0 
or, expressed in a non-dimensional form: 


L (—2B) r 
2 r —m 

Owing to the absence of exhaustive data relating 
to the variations of the value of the exponent and 
to the absence of any evidence regarding the variation 
of the value of the coefficient K within the range of 
the Reynold’s Number for airship forms in discussion, 
we are compelled to temporarily make use of the cor- 
responding values obtained for circular tubes. The 
results of tests obtained by Nikuradze' furnish us with 
an exhaustive idea of the variation of the value of the 
exponent and coefficient K within a large range of 
Reynold’s Number for circular tubes. The exponen- 
tial law expressing the velocity distribution in a tur- 
bulent boundary layer, as is well known, represents 
merely an approximation to the universal logarithmic 
velocity distribution. At the same time both the expo- 
nent » and the coefficient K decrease with the increase 
of the Reynold’s Number. Fig. 2 shows us the varia- 
tion of the values of 1/n and K versus the logarithm 
of the Reynold’s Number based on the length of the 
hull of the airship. 

It is clear that a diagram obtained in this manner can 
be used for calculation in the first approximation only. 

Further it will be shown that calculations based on 
Fig. 2 have been indirectly confirmed. 

Taking a separate case, we assume 


n = f (Re) = constant; K = 0.045 


and so on and then we obtain Clark Millikan’s equation* 
which at small Reynold’s Number gives good results. 


0.06 


/0 


8 Re 


Fic. 2. 


‘ Nikuradze, Gesetzmassigkeiten der turbulenten Strémung in 
glatten Roheren, Vorschungsheft No. 256, 1932. 


* Clark B. Millikan, The Boundary Layer and Skin Friction 


for a Figure of Revolution. Transactions of the A.S.M.E,, 


January 30, 1932. 
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Fig. 3 gives us an idea of the variations in the rela- 
tive thickness of the boundary layer in case of an 
increase of the Reynold’s Number. 

The thick lines represent the curves of the relative 
thickness of the boundary layer of the airship “Akron” 
hull calculated according to Eq. (6). 

Assuming n =f (Re) = constant according to Eq. 
(6) we shall obtain: 


+ = ¥ (Re) = constant (Re) "/*™ (7) 


For purposes of comparison the curves of the rela- 
tive boundary layer thickness computed from Eq. (7) 
for the same Reynold’s Number assuming » = 1/f are 
also given on Fig. 3. 

These curves are represented by thin lines. 

When comparing the relative thickness computed 
from Eqs..(6) and (7) for the same Reynold’s Num- 
ber we arrive at the conclusion that if the change of 
the coefficient K and the exponent are not taken into 
consideration the decrease of the relative boundary 
layer thickness becomes greatly exaggerated as the 
Reynold’s Number increases. 

When proceeding from model to full scale the abso- 
lute thickness of the boundary layer % will increase, 
owing to the increase of the dimensions of the body, 
whereas the relative thickness 8/L will decrease. 

When passing from model to full scale the part of 
the empennage operating in the boundary layer of the 
hull will decrease. 


% 
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This shows the greater efficiency of the full-scale 
tail surfaces as compared with the model ones. For 
the correct estimation of the forces and moments act- 
ing on the empennage in full scale it is therefore 
necessary to correct the normal forces and the pitching 
moment coefficients obtained when testing the empen- 
nage, in correspondence with the changes in the rela- 
tive boundary layer thickness. 

For the first approximation we can recommend the 
following method of introducing the necessary 
corrections. 

Assuming the velocity outside the boundary layer 
in the region of the empennage to be constant and 
equal to the mean velocity of the outside boundary of 
the layer (w, mean) we may calculate in what extent 
the efficiency of the part of the empennage lying out- 
side the boundary layer will be greater than that of 
the empennage lying inside the boundary layer. Let 
us evaluate the relation of the efficiency of the empen- 
nage working inside the boundary layer to that of the 
empennage outside the boundary layer by the coeffi- 
cient : 


ean)? d 


(us mean)? (us mean)? ia 2n+1 


Let us denote the area of the empennage operating 
outside the boundary layer by S, and the area of the 
empennage working inside the boundary layer by S, 
and introduce the notation of the effective area of 
empennage as: 

Sith S, (9) 


representing the area of some equally-effective empen- 
nage, working in a flow with a constant velocity along 
the section of the boundary layer. Knowing the thick- 
ness of the boundary layer in the region of the empen- 
nage for model and for full scale the values of S moaei 
and Sjurt scale May be computed from Eqs. (8) and 


(9) in model scale, i.e. 1/m’ Srsts scale Where m— is 
the linear transition coefficient from model to full scale. 
It is evident that 


1/m* Spuit scate > Smodet 


Hence it follows that in order to make use of the 
coefficients of normal forces and longitudinal moments 
acting on the empennage and determined by experi- 
menting with the model in full scale we must multiply 
the coefficient by the following ratio 


Sputt 


The thickness of the boundary layer 3 being known, 
the local friction force may be determined from Eq. 
(4). In case of a given magnitude of the frictional 
force +, the longitudinal aerodynamical load resulting 
from friction may be found. 

Practically, the part of the drag resulting from fric- 
tion can be considered as the total airship hull’s drag: 


THE BOUNDARY LAYER AND DRAG 


REVOLUTION 19 


OF A BODY OF 


L p y2 L ™ 
D -f = 2xK —— (=) F° 76"dz (10) 
0 2 v 0 


Reducing the drag to the non-dimensional volume 
coefficient C, and to the non-dimensional friction 
coefficient C, we obtain: 


Cn= - D = 9 K R 
L ( * ( L ) (11) 


where Vol. is the volume of the hull and 


Cy = — __ ["F ( a(+) 


where f is the surface of the hull. 
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Fig. 4 shows the variation of the logarithm of the 
drag coefficient Cp computed from Eq. (11) versus 
the logarithm of the Reynold’s Number for the 
“Akron’s” hull. For purpose of comparison the 
dependence of the logarithm of the drag coefficient 
Cp on the logarithm of the Reynold’s Number, assum- 
ing n= f (Re) =constant=1/F is also shown. As 
is well known this dependence is represented by a 
straight line. The end of the test curve determined 
by means of testing the model of the “Akron’s hull 
in the American variable density wind tunnel* is like- 
wise plotted on the diagram. 

We see that at full-scale Reynold’s Number (which 
interest us most) the coefficient C, determined under 
the assumption that f (Re) = constant = 1/f, 
i. e. according to Millikan’s theory greatly differs from 
the coefficients determined from Eq. (11) according to 
the theory in question. 

3 Hugh B. Freeman, Force Measurements on a 1/40-Scale 


Model of the U. S. Airship “Akron”, N.A.C.A., Report No. 
432, 1932. 
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The most convincing confirmation of the possibility 
of realizing the suggested theory in practice will be 
given by comparison of the friction coefficient C, deter- 
mined according to Eq. (12) and the friction coeffi- 
cients for flat plates determined according to Karman’s 
and Schlichting’s latest formulae deduced on the basis 
of the universal logarithmic law of velocity distribu- 
tion in the boundary layer. 

The results of this comparison are shown on Fig. 5. 

The curve C, for the “Akron’s” body computed from 
Eq. (12) bears the same character as the curves for 
the flat plate and is located somewhat higher. 

It can be assumed that at smaller fineness ratio 
(length to diameter ratio) the curve of friction 


coefficients C, will lie somewhat higher owing to an 


bts obtained by Hemp} fora flot 
© plate placed in 


of theoretical friction intensit 
distubution the mocke at 


D452 
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increase of the curvature and of the value of the static 
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SUMMARY OF RESULTS 


OR stiffeners or struts built up of flat plate elements 

and subjected to axial compression, theoretical calcu- 
lations can be made with a degree of accuracy sufficient 
for practical application, and will be on the conserva- 
tive side. Failure in such a member will occur when 
the applied stress is equal to the critical stress value 
in Euler column failure, torsional failure, or flat plate 
buckling failure, which ever is the smallest. The 
formulae and diagrams obtained (Fig. 5) enable the 
designer to determine the type of failure occurring in 
a practical case. 


STATEMENT OF PROBLEM 


Struts built up of flat plate elements are extensively 
used as stiffeners in aeronautical construction, and con- 
sequently knowledge of their strength under compres- 
sive load is of importance to the industry. All shapes 
of stiffener sections made up of thin flat sheet of rea- 
sonable length, may be subject to three classes of 
failure, namely : : 

(1) Local wrinkling or buckling of the sheet in the 
component parts, 

(2) “Buckling”, i.e. buckling of the stiffener as a 
whole acting as a long slender column in compression, 

(3) “Torsion Failure” due to torsional instability 
under compressional load. 

The question at once arises as to the points of transi- 
tion from one type of failure to another, i.e. the 
limits within which each type of failure may be expected 
to occur. It is apparent that any particular stiffener 
will experience that type of failure for which the 
critical failing stress has the lowest value. 

An investigation of angle section stiffeners of equal 


leg lengths has previously been made by A. B. Vossel- 
ler and C. C. Jerome,’ consequently as an extension 
of their research, the extensively used channel section 
was chosen for this investigation. 


THEORETICAL CALCULATIONS 


In order to present the results of the investigation 
in a clear and concise manner, it is necessary to express 
the formulas for the three different types of failure in 
terms of the In addition, it is 
desirable to have these parameters dimensionless. The 
formula for the critical stress causing buckling of a 
flat plate under compression was selected as the typical 
equation. This equation has been given in a recent 
paper by Donnell? in the form: 


same parameters. 


Ger = KE (thickness /width)? (1) 


where K is a dimensioniess coefficient depending upon 
the ratio of the length to the width of the plate and 
upon the edge conditions. The edge conditions in 
the problem under investigation result in K becoming 
a function of b, w, L, t and p» 
where == length of leg (depth of channel) 

w == width of back 

L = length of channel 

t == thickness 

== Poissons ratio (Taken as 0.3 for material 
tested ) 

as illustrated in Fig. 1. 


1A. B. Vosseller and C. C. Jerome, /nvestigation of Relation 
between Euler and Flat Plate Buckling of “L” Section Struts, 
Transactions of the A.S.M.E., Applied Mechanics Division 
August, 1934. 

2L. H. Donnell, The Problem of Elastic Stability, Transac- 
tions of the A.S.M.E., Vol. 5, No. 4, p. 4, 1933. 


Fic. 1. Channel 
Dimensions. 


Fic. 2. End Blocks Showing Grooves for Channel Ends and Hemispherical Ball. 
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(a) Euler Failure 

The theory for the failure of struts as Euler columns 
is well known and has been checked closely by von 
Karman,’ who originated the general method of vary- 
ing end conditions to counteract initial eccentricities 
which was used in this investigation and that of Vos- 
seller and Jerome.!| The conditions under which the 
stiffeners were tested were made as nearly pin ended 
as possible by mounting the channels in grooves on 
steel plates which rested on a hemisphere of hardened 
steel (cf. Fig. 2), hence the end fixity was taken as 
unity for all calculations. The effective length of the 
channel to be used in the Euler column formula, how- 
ever, exceeds the actual length due to the addition 
of the supporting blocks between the ends of the strut 
and the hinged end or point contact. Von Karman 
deduced that the exact value of the P gure, for such a 
strut may be computed from 


cot VP/EI (L/2) = a ¥P/EI 


where a = the additional length at one end due to the 
rigid supporting block. 

In this case, this leads to the result that the effective 
length is the original length plus the depth of the sup- 
porting blocks at both ends.* A standard channel length 
of 22 inches was used in this investigation, and the 
combined depth of the two end blocks from the bot- 
tom of the grooves to the top of the hemispheres was 
determined as 2.3” giving an effective Euler column 
length of 24.3” to be used in the Euler formula. 

In order to determine the value of A for the Euler 
column formula, the Euler critical stress is equated 
to the critical stress of the type Eq. (1) which gives 

E 7? (r/L)? = EK (t/b)? 


where y = radius of gyration of the stiffener channel 
section. Since the channels are of thin sheet, and the 
area is approximately (2)+ w)t, a simple calcula- 
tion gives b’/(2b + w) for the distance of the C.G. 
from the back, when ¢ terms are neglected. Calcula- 
tion of the radius of gyration r and substitution in the 
above formula gives 

Kewer = [(b/w)* + 2 (b/w)] 

/ {3[2(b/w) + 1): 


about the axis parallel to the back and 


Keuter = 7? 0* [1 + 6 (b/w) ] 


{12 (b/w)? [2 (b/w) + 1] } (2-b) 


about the axis perpendicular to the back. 


(b) Plate Failure 
The critical stresses for the flat plate failure of the 


3Th. von Karman, Mitteilungen tiber Forschungsarbeiten, 
Ver. Deuts. Ing., Berlin, 1909, 

4 Proof of this may be found in the Master's Thesis of the 
authors at the California Institute of Technology. 
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legs and back must be considered separately. The 
formula for the legs is given by an equation of the 


type 
Scr = KE (t/b)? 


in which b= length of leg. The critical stress for 
flat plate failure of the back is expressed by the same 
formula with width of back “w” substituted for “b”. 
If the leg is considered as having simple support along 
the edge where it joins the back, the value of Donnell’s 
K for this condition? is 0.4. Considering the same edge 
conditions for the back, which would result in simple 
support on all 1our sides, we obtain K = 3.6. When 
the edge joining the leg and back is considered as fixed, 
K becomes 1.2 for leg plate failure and 6.3 for back 
plate failure. (For the purpose of comparison, the 
value of AK used by Donnell in Ref. 2 is related to the 
K found in Timoshenko® and other books on elasticity 
by the equation Kp= KyX (1—p*) where 
is taken as 0.3). It is apparent that the exact value 
of AK for the specimens tested should be somewhere 
between these two extremes of simple support and fixed 
edge conditions, and must be determined by combining 
the effect of the legs and back on each other. 

The failure of the legs as plates buckling under com- 
pressive load is essentially a problem of the elastic sup- 
port given the legs by the back as the width of the 
back varies. It is, however, still a problem of plate 
buckling and may be analyzed in one of two or three 
methods all quite complex. Of these the Energy 
Method was taken as the least involved, and the fol- 
lowing theoretical expressions were developed by von 
Karman. 

The internal energy stored up in the plate during 
deformation is first set up using sine curve approxima- 
tions to the true deflection curves in the length direction 
and in the direction of the width of the back, assuming 
that the elements of the leg length b remain straight. 
The work done by the external force is then set up in 
terms of the deflection, and considering that the energy 
change is zero, the work done is put equal to the energy 
stored up in the plate. The critical value of the com- 
pressive stress for buckling is so determined and the 
value of the coefficient K for this type of failure found 
in terms of the parameters of the channel. 

Considering that the channel is made up of three 
thin plates simply supported or hinged at their top 
and bottom edges and deformed by the compression 
load, the energy of bending and twisting for the 
deformed channel is (cf. Fig. 1) 


Et L 1 &z\2 
&z 
where z is the deflection normal to any of the plate 


elements. It is assumed that the deflection of the leg 


5 S$. Timoshenko, Strength of Materials, Part II, p. 605, 
D. van Nostrand & Co., 1930. 
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may be represented by 
Z, = Av, sin NT y/L 
and that of the back by 
Z, = B (sintz,/w) (sin Nry/L) 
The angle between the back and leg is assumed to 
be preserved so that at the joint @2:/@x,= 22./81r: 
Then the total internal elastic energy is 


The external work done during buckling is given by 


w = f f ” (62/by)* dedy 


Equating these two energies we obtain ¢, and intro- 
ducing the condition that ¢ be a minimum, we get ¢-,,, 
Referring to Eq. (1) the result may be put in the form 
K =2n?V (b/w)? [((V +1) + 4V (b/w) 

x (1 — p) 1/12 — p?) (3) 


where V= 3/4 (b/w)* + 3 

It is realized that the assumption that the leg remains 
straight in the length direction is somewhat in error 
and will lead to too large values of K. However, to 
approach the actual condition any closer, an analytical 
expression must be chosen for the leg deformation, 
and the calculation becomes too involved to be pre- 
sented here. 


Torsion Failure 

The theory for failure of this type of stiffener due 
to torsional instability has been treated recently by 
H. Wagner and W. Pretschner.6 This method was 
chosen by the authors as being most suitable for their 
use and expresses the value for the critical load of a 
section under compression as 

P = Glo + Cpr/L*)/i's, 


where is, = polar radius of gyration of the section 
with respect to the shear point or elastic 
centrum, 
G@ = modulus of elasticity in shear, 
E = Young’s modulus of elasticity, 
Ie = effective moment of inertia for thin sheet open 
sections 
Cpr =constant determined by the bending rigidity 
of the flanges. 
This formula is made up of two parts which will be 
treated separately and combined later. 
or = P/A = [E [2/6 (1 + 


from which 
K, = 0/6 (1 + p) tp (4) 


®*H. Wagner and W. Pretschner, Verdrehung and Knickung 
von Profilen, Luftfahrtforschung, Band 11, No. 6. 
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Fic. 4+. Area Swept by 
Radius Vector for Cal- 
culation of 


‘ 


Fic, 3. Shear Distribution. 


In the determination of i°,, the shear center or elastic 

centrum is located’ at (cf. Fig. 3) 
e = W? bt/4 T,, 

then 2 
14+ 6(b/w) b/w[2+ (b/w)] 
12(6/w)*{2 (b/w) +1] * 3[26/w)+1P 

_ 


For Part II 

P, = E C,,/L? ity 
in which Wagner defines the quantity Cgr by the 
integral 

C = w? dA 


and shows that # = twice the area swept through by 
a radius vector from the shear center to the outline of 
the section. The physical meaning of w is given by 
the fact that in the case of torsion without axial 
stresses the axial displacement of an arbitrary point 
of the cross-section is proportional to w, hence the 
w distribution determines the warping of the cross- 


section. 
For the back 


For half the back then 


f wdA =1 = te-w*/24 
and for one leg 


=t f [(ew/2) — (wx/2)P dz = 
o 
[tw?/4] [eb — e b? + b°/3) 
Then for the total section 
Car = [2¢w?/4] [e?w/6 4 eb — eb? + 0/3] 
Substituting for e and putting this value of Cgrin the 
equation for P, gives a value which from o = P,/A 
gives the value of 
Ky = [2/12] (b*/L? [w/b] 
X [3 (b/w) + 2]/[12 (b/w)? +8 (b/w) +1] (5) 


and the total K for torsion = K,+ K, 


7S. Timoshenko, Strength of Materials, Part I, pp. 193-194, 
D. van Nostrand & Co., 1930. 
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It is seen from the formulas for K pyser and K qorsion 
that K is a function of b/w and of the parameter 
b*/Lt. Values of the latter corresponding to the differ- 
ent specimens investigated are indicated in the figures 
where K is plotted as a function of b/w using the 
above equations. A pige being merely a function of 
b/w is not affected by the value of 6*/Lt and is the 
same in all the figures. 


EXPERIMENTAL INVESTIGATION 


The determination of the compressive strength of 
channel sections involves three variables, namely thick- 
ness of material, width of back, and length of leg. In 
oider to determine the effect of each variable, it 
becomes necessary to vary only one at a time and this 
was done by using two thicknesses of material .025’ 
and .051” widths of back from 0.4” to 2.0”, lengths 
of leg from 0.4” to 1.5”. The tests were made with 
equal length channel section struts 22” long which, 
with the supporting blocks, etc., gave an effective length 
of 24.3” for Euler failure. 

Since accurate knowledge of the end fixity is essen- 
tial for the determination of column buckling loads, 
the apparatus used was the same as that described by 
Vosseller and Jerome’ except that the end plates were 
machined with perpendicular grooves in order to receive 
the channel specimens (cf. Fig. 2). Hardened steel 
discs were made for the ball to rest on and the discs 
shifted after each test. It was found that the slight 
support given a channel by the Ames dial gauge used 
for lining up the specimen would prevent Euler failure 
at the theoretical point and carry the load, in many 
cases, far above the Euler value. For this reason, 
after the channel was lined up (initial eccentricities 
removed), the dial gauge was removed before con- 
tinuing the test. 

The specimens tested were of 17ST aluminum alloy 
22” in length. Extreme care was necessary in squar- 
ing the ends, due to the tendency to have torsional 
failure from uneven stress distribution if the entire 
end of the channel was not flush in the bottom of the 
grooves. This was accomplished by milling off 
approximately .001” on each end after the ends had 
been chamfered. The channels were made with the 
grain (i.e. direction of rolling) of the material paral- 
lel to the angles. This was done due to the marked 
difference in the stress strain curve taken normal and 
parallel to the direction of rolling found in tests con- 
ducted at the Bureau of Standards for the N.A.C.A. 
on duraluminum sheet,* and to the common practice of 
manufacturing stiffeners with the grain in this direction. 

In order to properly line up the specimen in the 
cage, wooden jigs were constructed which would place 
the two hemispheres directly over one another when 


SE. E. Lundquist, The Compressive Strength of Duraluminum 
Columns of Equal Angle Section, N.A.C.A. Technical Note 
No. 413, 1932. 
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the end blocks were fitted in the jigs. An Ames dial 
gauge was then placed at the center of the specimen 
and the end plates adjusted equally, as load was applied 
until all the eccentricity of the channel was removed 
as indicated by a constant reading on the gauge for 
a considerable part of the expected load. The removal 
of all initial eccentricity must be very exactly carried 
out for satisfactory results due to its influence both 
in Euler and torsion failure. The gauge was then 
removed and the load increased by small increments 
until the specimen failed. 


DISCUSSION OF EXPERIMENTAL RESULTS 


The first tests made were of the ¢=.025 inch 
b = 04 inch channels varying the back from 0.4 inch 
to 2.0 inch, for which quite consistent Euler failures 
were obtained and which followed the Euler curve 
very closely (Fig. 5a) showing that the end conditions 
were as desired (fixity approximately unity). The 
two Euler curves correspond to failure parallel to the 
back and perpendicular to the back. The next set of 
channels tested was t = .025”, b = 0.6” which gave 
three torsion failures with the shorter backs and the 
remainder Euler failure as may be seen on Fig. 5b. 
It will be noted that the torsion failures obtained dur- 
ing the investigation were all consistently higher than 
the corresponding theoretical curves. The theoretical 
curves are based on an assumption of free ends, i.e. 
the ends are allowed to warp when the channel twists. 
This condition naturally is practically impossible to 
obtain and resulted in the channels resisting the tor- 
sional moments until a higher load than the theoretical 
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Fic. 5. Theoretical Curves and Experimental 
Test Points. 
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value was obtained. The channel would then fail in 


torsion and in failing would cause a local buckling of 
one corner of a leg indicating the tendency for the 
legs to warp in torsion failure. If, however, the stress 
reached the Euler value before this occurred, Euler 
failure would result, as clearly shown by this set of 
points following the Euler curve after three torsional 


failures. 

It being desirable to obtain plate failure also, chan- 
nels were chosen ¢ = .025”, b = 0.8”; t = .025”, 
1.5”: to insure plate failure of the legs. The 
test values obtained for these channels, as may be 
seen in Figs. 5c to 5f, followed the general curva- 
ture of the theoretical plate failure curve except for 
those points where b/w-> 1. It is suggested for design 
that the dotted curve faired through the experiment 
points be used. 

In this region, considerable scatter is apparent due 
to the tendency toward torsional failure, five of the 
failures being definitely torsional. It is felt that sev- 
eral of the other points were probably a combination of 
torsion and plate failures, since a small additional load 
after plate failure in this region would result in a 
twisting of the channel, but the initial failure could 
not be definitely designated as other than plate failure. 
For all the points shown, not marked as torsional, in 
this group of channels plate failure of the legs occurred. 
The low values obtained for this set of channels as 
compared to the theoretical plate failure curve is ex- 
plainable by the fact that in the theoretical analysis for 
this case, the assumption was made that the energy of 
the deformed section could be obtained by assuming 
that the legs deform as a straight line. This assumption 
will hold for very large values of zw, or small b/w ratios, 
but as w becomes small, or b large (large b/w ratios), 
then the assumption of a straight line deflection for the 
legs breaks down. This must be true since at the limit 
(w > 0, b/zv-> ©), the legs of the channel will act as 
plates fixed along one side, and will obviously not have 
a straight line deflection curve. Thus for large values 
of b/w the treatment previously given leads to a higher 
value of K than would be expected from experimental 
tests. Probable discrepancies in the end conditions 
may also be held partially accountable for the low 
value. 


CONCLUSIONS 


It is fully appreciated that the method of securing a 
channel strut in a structure will vary from no support 
to rigid support for the back, and that the degree of 
this support will greatly influence the amount of fixity 


given the edge of the legs. It is hoped, however, that 


this investigation will at least give some basis for 
making a first approximation of the strength of chan- 
nel section struts. 

The authors believe that 
demonstrated the feasibility of analyzing all sections 
built up of flat sheet elements, theoretically with some 
degree of accuracy, fundamentally they all 
reduce to the same question of determination of K. 
It will be noted that the experimental points check the 
theoretical values for the plate failure curve quite 
closely in the usual range of channel 
i.e. 0.4 < b/w < 0.7 and beyond this range the torsion 
failure region is entered and failure occurs in this 


this investigation has 


since 


sections, 


manner. 

Probably the most immediate benefit to be derived 
from this investigation is the derivation of values for 
K for all three types of failure which are expressed 
in terms of the dimensions of the channel. By compu- 
tation of the three A values by the formulae given 
and the selection of the minimum value, the engineer 
can quickly determine the type of failure to expect 
and can then compute his allowable stress for that 
type of failure. 


Fic. 6. Testing Machine Showing Typical Plate 
Failure of a Specimen. 
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Longitudinal Potential Flow about an Arbitrary Body of Revolution with Application 
to the Airship “Akron” 


Ricuarp H. Situ, Massachusetts Institute of Technology 
(Received April 1, 1935) 


SUMMARY 


The flow potentials for axial flow about any arbitrary 
surface of revolution of easy form are found by a new 
use of the classical ellipsoidal harmonic solutions of La- 
place’s Equation in three dimensions. From the flow 
potentials the surface velocity and pressure are computed 
in the usual way. The treatment is applicable to the solu- 
tion of the potential flow about a number of bodies of 
practical importance. For example, it is equally applicable 
to airships, to fuselages which are bodies of revolution, 
to engine nacelles or to submarines which are deeply sub- 
merged. 

In order to compare the theoretical and experimental 
pressures in a specific case of interest, the surface pres- 
sures over the U. S. Airship “Akron” are computed and 
compared with the Langley Field measurements on a 1/40 
scale model. The agreement is excellent. The longitudinal 
inertia coefficient is also found, in the case of the applica- 
tion to the “Akron,” from the flow velocity potential and 
added to the effective coefficient of the boundary layer 
as determined from the Langley Field results. The total 
coefficient so found agrees well with values previously 
assumed. 


INTRODUCTION 


OTABLE progress has been made during the 
+ last few years in solving one of the difficult 
problems of classical hydrodynamics, namely, the prob- 
lem of: finding the potential flow about arbitrary bodies 
of easy form. The first major advance in this field 
was made by von Karman,' who found the approximate 
flow velocity and pressure about an arbitrary airship 
hull, for both axial and inclined motion in ideal air. 
The principle underlying the von Karman method is due 
chiefly to Rankin who was the first to substitute for 
the surface of a solid body the surface of separation 
which divides a source to sink or a doublet flow from 
the flow of the outer uniform stream when the two are 
superimposed. In the von Karman application of this 
principle to the arbitrary airship hull, two separate 
solutions are made and combined, namely, the solution 
for axial flow in which the hull is replaced by the 
equivalent system of sources and sinks properly dis- 
tributed along the hull axis and the solution for normal 
flow in which the hull is replaced by the equivalent 
system of doublets similarly distributed. Both solu- 
tions are in turn combinations of two separate solutions, 
each for a half body of revolution. The crux of the 
method consists of the determination of the several 


1 Th. von Karman, Berechnung der Druckverteilung an Luft- 
schiffkarpern, in Abhandlungen aus dem Aerodynamischen In- 
stitut on der Technischen Hochschule, Aachen, Germany, 1927. 
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unknown strengths of the source and sink or of the 
doublet segments, as the case may be, into which the 
half body axis is arbitrarily divided, so that the surface 
of separation will coincide with the surface of the 
arbitrary half body. The strengths are determined by 
an ingenious solution of a system of n linear equations 
where 1 is the number of source or sink or doublet seg- 
ments expressing the boundary condition, Y = 0, in 
terms of the unknown strengths and the coordinates 
of the surface. When the strengths of the several 
source or sink or doublet segments along the axis are 
known, the velocity potential of the flow is known and 
the velocity and pressure over the surface of the half 
body can be easily determined. In this way von Karman 
developed first the flow potential and pressure about 
two half bodies of revolution, one corresponding to 
the bow half of the arbitrary airship hull and the other 
to the stern half. He then derived the hull form and 
the curve of pressure distribution over it by cutting 
each half body and its pressure curve to the correct 
length and joining them amidships. Certain approxi- 
mations and assumptions incident to the splicing were 
made and investigated. 

The analogous development in two dimensions was 
made by the author who found the axial flow potentials 
about an arbitrary strut by a similar use of a system 
of parallel line sources and sinks lying in the strut 
plane of symmetry. In this study the arbitrary body 
was treated as a single closed surface, not as two 
separate half-bodies, and the whole system of sources 
and sinks was determined by a single solution of a 
system of thirteen linear equations.* 

The next major advance in this field consisted of a 
series of studies on the two-dimensional potential flow 
about arbitrary wings. These studies which began with 
two or three treatments of the special cases of thin 
sections, have now been sufficiently advanced through 
the recent treatment by Theodorsen to include a satis- 
factory analysis of the potential flow about wing sec- 
tions of quite arbitrary form.* Since these wing 
studies, however, depend chiefly upon the method of 
conformal transformation of elementary flows whose 
potentials are already known, they are of no direct 
interest here. 

The present study returns to the original problem in 
an attempt to solve more satisfactorily the potential 


2R. H. Smith, Aerodynamic Theory and Test of Strut Forms, 
N.A.C.A. Report No. 311, Part 1, 1929. 

3Th. Theodorsen, Theory of IWing Sections of Arbitrary 
Shape, N.A.C.A. Report No. 411, 1932. 
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flow about a body of revolution of easy but otherwise 
arbitrary form. The treatment first develops the axial 
flow potentials about any such body directly by a new 
use of the classical zonal harmonic solutions of La- 
place’s equation in three dimensions which are already 
known. From these potentials, the velocity and pres- 
sure distribution and the ideal virtual mass of the 
arbitrary body are then determined. Finally, in order to 
illustrate the method, a complete application is made 
to the U. S. Airship “Akron”, although the method 
is equally applicable to airplane fuselages and engine 
nacelles, to submarines and to other bodies of practical 
usefulness whose general forms are surfaces of 
revolution.° 
THEORETICAL TREATMENT 

The method of the present investigation is based 
upon the zonal harmonic expressions for the velocity 
potential and the stream function for axially symmetri- 
cal flow, as given by Lamb,® 


= Pn X Qn 


k 


which are given in terms of the spheroidal nondimen- 
sional coordinates ; 


where b = 1 — + /k’]. 
system, all meridional planes are divided alike by 


d Qn (5) 


In this coordinate 


orthogonal families of confocal spheroids, € = const., 
and hyperboloids, » — const., whose common focii are 
atx = +k. The relation between the two potentials 
is then, 

r O8p r O8¢ 


where r is the radial distance from the axis of sym- 
metry and where 6s, and és, are elements of distance 
normal to the local spheroid and hyperboloid respec- 
tively. The P, and Q, functions are the first and second 


kind of zonal harmonics of the n™ order. Since each 


re or value of 7 in Eq. (1) gives a separate solution 
f Laplace’s equation it follows that the sum of a 


“In rm application to the “Akron” Airship an approximate 
determination of the viscous part of the virtual mass is also 
given in order to be complete, although it hardly belongs in a 
study of potential flow. 

‘A similar theoretical treatment of the axial flow about an 
sdhdimaase body of revolution, but without application to an actual 
airship, has been developed independently by Carl Kaplan, at 
the Langley Memorial Laboratory. His studies may be found 


in N.A.C.A. report No. 516 on “Potential Flow about Elongated 
Bodies of Rez volution,” 

being printed. 

mae” Lamb, Hydrodynamics, 5th ed., pp. 131-132, Cambridge, 
IG 


which appeared while this paper was 


series of solutions obtained by varying the order of 
P, and Q, will also give the potentials of a possible 
flow. The potentials of such a flow are, 


= (u) Qn 


i 


y'= +1) G— — ——— 


1 


Furthermore ¢' and y' in Eqs. (3) can be multiplied 
as usual by any arbitrary constant or increased or 
decreased any constant amount without affecting their 
validity as potentials of a possible flow. Use will be 
made of these constants to meet the necessary boundary 
conditions at the surface of the arbitrary body of 
revolution.* 
The potentials of a uniform stream flowing with the 

velocity, + U, are 

oo = Uxr= 

vo = = Uk? (1 — — 1)/2 (4) 
If a multiply factor UkA,/2 is inserted in Eqs. (3) 
and if the constants in Eqs. (4) are added, one obtains, 


1 


An  dPniu) dQn 


1 


(5) then give the potentials for the axial flow 
about a body of revolution which is at rest in a fluid 
moving at infinity with the uniform velocity + U. 

The unknown coefficients, A,, in Eqs. (5) will now 
be determined from the condition that the stream sur- 
face, ¥ = O shall coincide with the surface of the 
arbitrary body of revolution. In order to meet this 


condition the equation, 


dQn(t) _ 
+) n(n + 1) (6) 


du dt 


1 


must hold at all surface points except at the bow and 
stern where 1 — p? = O*. The harmonic derivatives 
which appear in Eq. (6) are geometric terms which 
have a definite value at each point on the arbitrary 


surface. Before the harmonic derivatives can be 


evaluated, however, the hull must be suitably oriented 
with respect to the facii of the spheroidal coordinate 


*One notes in passing that the two potentials in Eqs. (1) 
and (3), and hence the corresponding flow velocity vanish at 
infinity since both QO, and its derivative are zero at infinite 
values of ¢. 

* Because of this restriction the subsequent treatment will of 
course, be invalid at these two points. This defect, however, is 
not serious since the pressures there are already known. 
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system. A similar argument to the one followed by 
Theodorsen in the development of the two-dimensional 
flow about arbitrary wings,® leads to the location of 
the focii at the points midway between the bow or 
stern of the hull and the respective centers of curva- 
ture of the hull surface at the bow and stern, Fig. 1. 
When so oriented the hull departure from the nearest 
spheroid is a minimum. 

After the hull focii have been determined k is known 
and the (y, ¢) coordinates of any number of points on 
the hull surface can be computed from Eqs. (2). One 
can then evaluate n orders of the products of the zonal 
harmonic derivatives in equation (6) at each of m sur- 
face points and so obtain a system of m linear equations 
which, when solved, yields the m unknown coefficients 
ri er A,. The first of Eqs. (5) with these 
values of A, then gives the velocity potential of the 
axial flow about the arbitrary hull. It is noteworthy 
that if the first order only is present the other orders 
vanishing, the velocity potential reduces to that given 
by Lamb and others for the axial flow about a simple 
prolate spheroid. 

Having now the velocity potential for axial flow about 
the arbitrary hull, the velocity and pressure distribution 
can be computed without great difficulty. If as before 
3s, and &s, designate elements of distance normal to 
the local hyperboloid and spheroid respectively, then 
after Lamb, 


The local resultant velocity may be written 


v= y+ (ir) 


ab OF 
Ou 


Since the two velocity potential derivatives can be found 
immediately from Eqs. (5), the velocity can be evalu- 
ated from Eq. (7) for any surface or field point (pn, ¢) 
desired. The corresponding pressure relative to the 
stagnation pressure, p,, follows easily from the equation, 
=1— (4/0)? (8) 


P/Po 

The part of the virtual mass, m,, of the arbitrary 

surface which is associated with axial potential flow, 

can now be derived from the usual equation of field 
10 


energy, 
27 = m,U?=—p for aa (9) 


where x and a are the surface normal and the element 
of surface area respectively. The integration of Eq. 


* Lie. 


Hydrodynamics, 5th ed., p. 118, Cambridge, 1924. 
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10068 


—Ue 


Fic. 2. 


(9), unfortunately, cannot be carried through analyti- 
cally without some difficulty. It can, however, be 
carried through graphically very simply. In order to 
evaluate 2T graphically, one first assumes for con- 
venience that the distant fluid is calm and that the body 


moves through it with the forward axial velocity — U. 
The velocity component normal to the surface will 
then be 0¢’’/0n = — U cos a, where a is the inclination 
of the surface normal to the axis of symmetry. Sub- 
stituting in Eq. (9), one obtains. 
m, U?=pU | cosada 
(10) 


=2npU f cosa ds 


the elements of surface area, da, being 2zr,ds. The area 
under the curve formed by plotting the surface value 
of $’’r,cos a against the corresponding surface posi- 
tion, s, where s varies from zero at the bow to the de- 
veloped length of the surface at the stern will then give, 
when multiplied by the proper factors, the virtual mass 
of the potential flow outside the boundary layer. 
The corresponding inertia, coefficient, K,, is obtained 
by dividing m, by the mass of the displaced fluid, 
namely, 


= m,/p Vol. (11) 


It is also possible to define a virtual mass for the 
flow within the boundary layer provided the energy 
of the boundary layer is proportional to the square 
of the surface speed. The latter condition is fulfilled 
if neither the local velocities at various levels in the 
boundary layer relative to the surface speed nor the 
layer thickness vary appreciably with the surface speed. 
Experimental measurements will be referred to pre- 
sently which indicate that this is approximately the 
case, at least for bodies of revolution of easy form 
such as an airship in axial motion. 
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The energy of the boundary layer in any case, is, 


STERN 6 
B= fas f (r+ y 008.0) (Us; —u) *dy 
BOW 0 


where the terms are as defined in Fig. 2. 
If the boundary layer energy is proportional to 
square of the surface speed, U,,? Eq. (12) yields, 


the 


STERN 


pVol. pUsVol. J Vol. 3\ Vol.2/3 
BOW 


6 
2 
x cos a) (:- dy (13) 
Us 


where the constant, K,,, is the longitudinal inertia 
coefficient of the boundary layer, and the constant, 
Mr, 18 its virtual mass, and where Vol. is the volume 
of the hull. Eg. (13) can be integrated easily by 
graphical methods if the form, size, and speed of the 
airship are given and if the local velocity distributions 
across the boundary layer at a sufficient number of 
surface positions are known. A difficulty arises, how- 
ever, in deciding how and where to sever the boundary 
layer from its wake. In the application that follows 
this difficulty is met by arbitrarily extending the curve 
formed by plotting the integrand in Eq. (13) against 
s to the stern by continuous extrapolation and dropping 
it abruptly to zero, Fig. 5. The error in the total 
mass which this extrapolation introduces is probably not 
appreciable. The area under the curve of Fig. 5 is, 
of course, the inertia coefficient for the boundary layer. 

The total effective mass of the airship in longitudinal 
motion is composed therefore of three parts,—the 
natural mass of the airship, the virtual mass of the 
exterior potential flow and the virtual mass of the 
boundary layer, namely, 


m = m,+m,+m, (14) 


where #1, is the natural mass. Similarly, the total longi- 
tudinal inertia coefficient, K, becomes 


K,=14+K, + (15) 


APPLICATION TO THE U. S. AtrrsHtp AKRON™ 


The treatment above which leads to the surface pres- 
sure distribution and the inertia coefficients of an 
arbitrary airship hull in axial motion will now be applied 
to the U. S. Airship “Akron,” Fig. 1. In order to carry 
out the application one must first lay out the hull 
outline as in Fig. 1 and locate the focii of the (p, ¢) 
system of coordinates at the midpoints of the two radii 
of surface curvature at the bow and stern. Since k 
is then known, being the half distance between the 
focii, the (nu, €) coordinates of the hull contour which 
correspond to the (x, r) coordinates can be computed 
by use of Eqs. 2. One needs next the expressions for 
the various orders of the two zonal harmonics and of 
their first derivates. Since these harmonics are not 
common they will be given here up to the fourth order.'* 


™ Those who wish to carry the application through figure 
by figure may borrow a copy of the tables from the Aero- 
nautical Engineering Division of the Massachusetts Institute 
of Technology. 

” Byerly, Spherical Harmonics and Fourier Series, p. 190. See 
also H. Lamb, Hydrodynamics, 5th ed., p. 140, Cambridge, 1924. 
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P,(p) =p 
Py (p) = 5 — 3 
P, (a) = 4 3/8 
Q (2) =34Q, Q, /2 
Q, (2) = (2)/38 —2Q, (2)/3 
9, (0) = 749%, (0)/4 — 3 9, (0)/4 
P’, (nu) = 
(p) = 1 
= 35 — 15 p/2 
=Q’, =-—I// 1) 
(0) = 3 40/2 —O/2 + 3 0,2 
(2) = 5 2 — 20/3 + 5 
= —3 O/4 + 7 Q,/4 


where, of course, 


| 


n = (d/dy)P,(p) and 
Q’, = (d/dZ)Q,, (2) 

These ten harmonics and their derivatives were next 
evaluated at sixteen surface points on the hull, Fig. 3. 
One can now choose the minimum number of sur- 
face points and their positions which appear to be 
necessary to determine the departure of the hull con- 
tour from the nearest spheroid. In the case of the 
“Akron” airship whose form departs only slightly 
from the nearest spheroid, or from the nearest curve 
of constant € in the spheroidal coordinate system, four 
surface points in addition to the bow and stern are 
quite sufficient. The four intermediate surface points 
chosen were f, hi, j, and /, Fig. 1. The harmonic deriva- 
tives at these points P’,, P’:, P’s, P’s,Q’1,Q’2, Q’s, and 
being already known, the following system of linear 
equations based on Eq. (6) can be set up in terms 

of the four unknown coefficients A,, 


1+ A, (P’, Q’,) A, (P’, Q’,) 
+A 


3 (P’, Q’;) + A, (P’, Q’,) 7/20 = 9 
] A, (Pe Q’,)1/2 A, Q’.),/6 
+ A, (P’,Q’,)n/12 + Ay ,/20 = 0 as) 


1 + A, Q,) ;/2 + A, (P’, Qe) ;/6 
+ A, (P’, Q’,) j/12 + A, Q,) ;/20 = 0 


1 4.4, (P’, Q’,) 1/2 + Ay (P’, Qe) 1/6 
+ 


A, (P’, 1/12 + 1, (P’, Q’,) ,/20 


The solution of this set of equations is, 


A, = + 0.0661; A, = — 0.0386; (17) 


1 
A, = + 0.0166; A, = — 0.0388 


In order to prove that Eq. (6) in four terms is suff- 
ciently accurate as the equation of the arbitrary hull, 
test solutions were made at the five intermediate posi- 
tions using the four known A,’s. At all five of the 
test positions tried Eq. (6) was satisfied well within 
engineering precision. 

It is now possible to find the general velocity poten- 
tial, ¢, of the axial flow about the arbitrary hull from 


RECHARD Hf. 


SMITH 


the first of Eqs. (5), where » goes from 1 to 4, and 
A, has the value given in Eqs. (17). In the present 
application to the “Akron,” ¢ was computed for the 
16 surface points, Fig. 1, by evaluating the four 
harmonic products (P, Q,) at each point, multiplying 
each by its proper A, from Eq. (17), and summing 
with 2 €. 

The potential gradients (2/kU)(0¢/d,) 
(2/kU)(2/2%) can now be readily obtained from 
the equation of the velocity potential, Eq. (5), by 
simple differentiation, and evaluated numerically at the 
16 surface positions on the airship hull, excepting again 
the bow and stern positions for which the solution does 
not apply. When these values are inserted in Eq. 
(7) and the right side of Eq. (7) evaluated, one ob- 
tains the values of the velocity squared at each surface 
position. The relative surface pressures follow immedi- 
ately from Eq. (8) and are given in Fig. 3 where they 
are plotted against ship length along with the experi- 
mental pressure measured by Freeman.'* The agree- 
ment between the theoretical and experimental pressures 
is excellent except, as usual, near the stern where 
agreement is not expected because of the separation of 
the real air flow from the tail and near the bow where 
an unexplained irregularity appears in the experimental 
pressures. 

The velocity potential of the ‘“Akron’’ will next be 
used in a determination of the longitudinal virtual 
mass and the corresponding inertia coefficient A,, 
of the potential air flow outside the boundary layer. 
According to Eq. (10) 


and 


my = cos « (18 
where ds is an element of length lying in the surface 
along a meridian and where ¢” is the surface velocity 
potential of the airship moving through otherwise still 
air. From Eqs. (4) and (5), ¢”’/Uk in Eq. (18) 1s. 


6” 


and has the values which are given in Fig. 3 at the 14 
surface points between the bow and stern. Since its 
component parts are now determined, the integrand 
(¢"/Uk) (r,cosa/k), in Eg. (18) can_ be easily 
evaluated and plotted against s/k for integration, Fig. 
4. The value of the integral, which is the area under 
the curve, is found to be 0.001110. The virtual mass 
of the axial flow outside the boundary layer of the 
“Akron” is from Eq. (18), 


mp = 0.00222 = p k* = 0.0534 p Vol. 


where k®/Vol. == 7.67, and the corresponding inertia 


® Hugh B. Freeman, Measurements of oat w in the Boundary 
Layer of a 1/40-Scale Model of the U. S. Airship “Akron,’ 
N.A.C.A. Report No. 430, 1932. 
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coefficient, Eq. (11), is * 


K,, = ™,/pVol. = 0.0534 (19) 


As previously stated, the condition which makes it 
possible to define an analogous virtual mass for the 
flow within the boundary layer is that the energy of the 
boundary layer flow shall be proportional to the square 
of the surface speed, and this requires in turn that the 
boundary layer thickness and the local relative velocity 
within the layer shall not vary with the surface speed. 
It will be seen from Table 1 of Freeman’s report?’ on 
the measurement of the local velocities in the boundary 
laver of a 1/40 scale model of the “Akron,” that at 
most of the surface positions only small plus and 
minus differences were found in the relative velocities, 
u/U,, at the two higher airspeeds, and that no differ- 
ences were found in the boundary layer thickness. It 
seems safe to assume, therefore, that Eq. (13) is 
approximately valid. On this assumption, the integrals 


6 
f + y cos ( 1— dy were graphically 
3 


“Tt is of interest to compare the axial inertia coefficient for 
the “Akron” with that for the nearest spheroid. The fineness 
ration of the nearest spheroid is 5.9 whose axial inertia coefh- 
cient is K,=0.046. See H. Lamb, Hydrodynamics, 5th ed., 


p. 146, Cambridge, 1924. 
Tot. 


evaluated at each of the 11 stations chosen by Free- 
man, Fig. 3 of his report, by use of the data given in 
his table I for 7, y and u/U;and after multiplying 
by (U%s/U,)?, Eq. (13), are plotted 
against S/Vol.’* for integration, Fig. 5. The area 
under the curve in Fig. 5 which gives the inertia coeffi- 
cient K,, of the boundary layer, by Eq. (13), is found 
to be 0.0187 if one makes the simple assumption of 
cleavage between the boundary layer and the wake 
streams which has been referred to above.'* 

One finds in summary that the effective longitudinal 
mass and the corresponding inertia coefficient, of the 
“Akron” airship hull are by Eqs. (14) and (15), 
respectively, 

m = p Vol. (1 + 0.0534 + 0.0187 
K, = 0.0721. 


= 1.0721 p Vol. 


The value of the inertia coefficient which has been 


assumed at least in one instance is'® A, = 0.08. 


™S. J. De France and C. P. Burgess, Speed and Decelera- 
tion Trials of U.S.S. Los Angeles, N.A.C.A. Report No. 
318, 1929. 

* As this paper was going to print, the author learned that 
a similar evaluation of the effective inertia coefficient of the 
boundary layer of the Airship “Akron” had been made previ- 
ously by Dr. Klemperer of the Goodyear Zeppelin Corporation. 
His results, however, have never been published. 
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Book Reviews 
By R. Paut Harrincton 


Aerodynamic Theory, Dr. WILLIAM FrepERICK DURAND, 
Editor-in-Chief; Vol. III, Julius Springer, Berlin, 1935; 354 
pages, 167 figures. 

The third volume of this scientific work by leading authorities 
in aerodynamics ably maintains the excellent standards set by 
its predecessors in offering four additional divisions, three of 
which cover phases of fluid mechanics and one which treats 
of experimental methods and wind tunnels. 

Professors C. Witoszynski and M. J. Thompson discuss the 
theory of single burbling in division F, wherein, first for primi- 
tive circular profiles and later for other transformed profiles, 
an additional discontinuous term is added to the basic classical 
complex potential whereby the type of flow around the profile 
is altered. The result is the creation of a layer of discontinuity 
behind the profile which has its origin at a single point on the 
profile surface from which arises the expression “single 
burbling”. At the same time the calculated values of lift are 
decreased and are more nearly in accordance with the experi- 
mental values. 

Expressions for the air forces, including profile drag (less 
friction), and for the moment are derived. The general case 
of two dimensional flow and also the influence of a finite span 
are also treated. The division is an interesting and important 
contribution to aerodynamic theory. 

In division G, by far the largest of the present volume, Pro- 
fessor L. Prandtl has presented the mechanics of viscous fluids, 
a field in which he has pioneered for years. In addition to 
his usual excellent presentation of the earlier development of 
the laws and conditions governing the flow of real fluids, many 
of the more recent aspects are discussed. These deal mainly 
with advances in turbulent flow theories in regard to the origin 
and growth of turbulence, mixing length, and Professor von Kar- 
man’s development of the logarithmic velocity distribution for 
turbulent boundary layers. 

The recent work of Nikuradse and others on the effect of 
roughness on the flow in tubes and on plates is an interesting 
section. Turbulent shearing stress and velocity characteristics 
in flows accelerated or retarded by suitable pressure distribu- 
tions are presented according to the theories and experiments 
of Gruschwitz and or Buri. Prandtl concludes the division with 
a discussion of the resistance of bodies in turbulent flows and 
experimental methods for the determination of such resistance. 

The mechanics of compressible fluids are set forth in 
division H by Professor G. I. Taylor and J. W. Maccoll 
in a clear and concise manner. The propagation of small 
plane disturbances and latter of large disturbances and shock 
waves is dealt with. The conception of the shock 
wave and conditions within such a wave are_ treated 
with the apparent result that a permanent disturbance can be 
propagated under the action of conductivity or viscosity or a 
combination of the two. The means for calculating the velocity 
distribution within a shock wave is developed. The comparison 
is drawn between the dissipation of mechanical energy into 
heat which occurs with such a wave to the dissipation behind 
a flat plate normal to a stream of fluid. 

Shock waves may not occur when the speed of sound is 
exceeded in the case of steady flow through channels and jets. 
Such cases and relevant experimental results are described in 
the second chapter of this division. The two dimensional flow 
at high speeds less than that of sound is illustrated for a cir- 
cular cylinder and airfoils showing that for the latter case the 
relationship between circulation and lift is unaltered by com- 


pressibility. Lastly, the cases of plane flow at sunersonic speeds 


past thin airfoils, past a wedge and a cone, and around sharp 
or curved surfaces, and pressures in a Pitot tube or at the nose 
of a projectile are developed along with expressions for the 
air forces and moment for the former, which check the experi- 
mental results of Stanton. 

In division I, the last of the present volume, Professor A. 
Toussaint gives a historical account of the beginnings of aero- 
dynamic research, together with brief descriptions and develop- 
ments of various types of wind tunnels. The economics and 
energy relations inherent in each type are discussed in an inter- 
esting manner. The influence of the dimensions and boundaries of 
air-streams upon lifting systems and the corrections necessary for 
converting the experimental results into those for an infinite 
stream are derived according to the theories of Villat, Glauert, 
Terazawa, Carafoli, and others. 

Eastman Jacobs, in part 2 of the final division, presents a 
discussion of experiments conducted mainly at the N.A.C.A. 
on scale effect of such bodies as the flat plate, airship models, 
spheres, and airfoils and gives some interesting deductions there- 
from. Further descriptive material concerning the variable 
density wind tunnel at Langley Field is set forth together with 
recent experimental methods for the investigation of aerody- 
namic phenomena at high speeds. Numerous references to 
aerodynamic installations in various countries are included 
in the bibliography wherein descriptive material may be found 
regarding the form and construction of such installations and 
of their instrumental equipment. 


Handbook of Aeronautics, Caprain J. LAURENCE PRITCHARD 
and SqguaprRon LEADER C. G. Burce, Editors; Vol. I, Sir Issac 
Pitman and Sons, Ltd., London, 1935; 720 pages, 300 figures, 
$7.50. Vol. II, ANpREw Swan; Sir Issac Pitman and Sons, 
Ltd., London, 1935; 400 pages, 192 figures, $4.50. 

The second edition of this handbook finds a great part of the 
first edition (1931) rewritten and considerably enlarged. The 
volumes deal with modern aeronautical engineering practice and 
airplane engines respectively, and provides for aeronautical 
engineers a handbook similar to those available in other branches 
of engineering. The publication is under the authority of the 
Council of the Royal Aeronautical Society and the contributors 
are those whose names are in high standing in the aeronautical 
field. 

A considerable part of Vol. I is given to technical aerody- 
namics and performance, and has a copious supply of tables 
and curves of aerodynamical data applicable to everyday prob- 
lems in airplane performance and design. The remainder of 
the volume presents, in sufficient detail, descriptions, explana- 
tions, and data on construction, materials, meteorology, instru- 
ments, wireless, aerial photography, and design. Numerous 
references to less frequently used information is included in 
the various bibliographies. 

The second volume covers the design and practice in regard 
to aircraft engines, and presents a very comprehensive discus- 
sion and compilation of data under the various divisions per- 
taining to performance, design, carburetion, cooling, ignition, 
lubrication, engine tests, propellers, and others. Many diagrams, 
tables, and charts offer a very complete supply of information 
to the engine designer or experimentalist. 

The Handbook, as a whole, is written in straightforward 
fashion and seems to have sorted out the important engineering 
data and facts from the unimportant. The divisions are clean 
cut and numerous repetitions are avoided, a fact which in 
itself makes for a worthwhile handbook. 
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Institute Notes 


BENEFACTORS 


Through gifts of One Thousand Dollars, Donald W. Douglas 
and E. L. Cord have joined the growing group of members 
whose donations have made possible the work of the Institute 
and the publication of the Journal. Their names have been 


added to the list of Benefactor Members of the Institute. 


VISITORS TO THE INSTITUTE 


The following foreign members and guests of the Institute 
have visited the Skyport during the Summer: Dr. Carl Wiesels- 
berger, Professor and Director of the Aachen Aerodynamisches 
Institut, Aachen, Germany; Messrs. Henry Knowler, Chief 
Engineer, John Lord, Managing Director, and Captain Stuart 
D. Scott, Technical Advisor and Chief Test Pilot, of Saunders- 
Roe Ltd., Cowes, England; Stanislaw Krzyczkowski, Technical 
Manager of the Polish Airlines “LOT”, and President of the 
Aeronautical Engineering Association in Warsaw, Poland; 
Dr. Ing. Francois Misztal, Engineer of the National Aviation 
Establishment, Warsaw, Poland; Wilhelm Challier, of the 
Technical Institute of Warsaw, Poland; J. Van der Heijden, 
Inspector of Aircraft, The Hague, Holland; Squadron Leader 
D. C. Mackenzie Hume, of the Royal Canadian Air Force; 
Dr. G. Kaumann, of the Deutsche Luft Hansa and of the Asso- 
ciation of German Aircraft Manufacturers, Berlin; Mr. W. S. 
Shackleton, Aeronautical Consulting Engineer, London, Eng- 
land; Messrs. A. C. Barlow, Chief Engineer, W. Broadbent, 
Works Manager and Director, and R. T. Youngman, Technical 
Director, of the Fairey Aviation Company Ltd., Middlesex, 
England; Mr. A. Hessell Tiltman, B.S., F.R.Ae.S., Managing 
Director of Airspeed Ltd., Portsmouth, England; Mr. L. F. 
Bouman, Manager for Special Operations of the K. L. M., The 
Hague, Holland; Dr. Vilhelm Moser, Oberingenieur and Bevoll- 
maechtigen den Telefunken-Ges. fuer Drahtlose Telegraphie 
m. b. H., Berlin; Dipl. Ing. Gerald Klein, Siemens-Apparate u. 
Machinen G. m. b. H., Berlin, Germany ; Flight Lieutenant C. S. 
Staniland, Chief Test Pilot of the Fairey Aviation Company of 
Middlesex, England; Mr. Richard J. Moffett, Chief Engineer 
of Canadian Vickers Ltd., Canada, and President of the Com- 
mercial Aircraft Operators Association, Canada. 


University or MICHIGAN SUMMER SCHOOL 


Professor R. V. Southwell, M. I. Ae. S., of Oxford Uni- 
versity, gave a course of lectures on Applied Theory of Elas- 
ticity at the University of Michigan Summer Schoo! and, with 
Professor Timoshenko, also a member of the Institute, con- 
ducted a joint Seminar on Applied Mechanics. While not 
primarily devoted to aeronautical structures Professor South- 
well introduced new methods of calculation of stresses in 
redundant space frame work which may be of importance to 
aircraft designers. 


INTERNATIONAL CONGRESS FOR APPLIED MECHANICS 


The American Committee for the International Congress 
for Applied Mechanics met at Ann Arbor at the time of the 
meeting of the A.S.M.E. Applied Mechanics Division in June. 
The Chairman, Dr. J. S. Ames, who was abroad, was repre- 
sented by Dr. George W. Lewis. Other members present were 
Professors Timoshenko, yon Karman, and Hunsaker. 

The Congress will be held in Cambridge, Massachusetts, in 
September 1938, with Harvard and the Massachusetts Insti- 
tute of Technology as joint hosts, and the Institute of the 
Aeronautical Sciences and the A.S.M.E. cooperating. 


THe AERONAUTICAL INDEX 


The Works Progress Administration which has now taken 
over the direction of all projects such as the Aeronautical 
Index has approved an appropriation of $34,650 to be expended 
in carrying on the work of expanding the work already under 
This grant will provide salaries for thirty-two employees 
The personnel has been classified 


way. 
for the period of one year. 
so as to furnish a well-balanced group of research workers, 
translators, typists and clerks. 

Over 120,000 cards are now in the 
are added each week. 

It is the aim of the Institute, not only to index all books 
relating to aeronautics, but through a subject index file have 
the main subjects with which they deal classified so that the 
more accessible to 


Index and about 5000 


contents of volumes will be made 
workers in aeronautics. 

The Centre de Documentation Aeronautique Internationale 
which publishes, bi-monthly, Aeronautiques, a 
classified digest and index of all aeronautical books, articles, 
papers, reports, etc., have sent a complete set of these valu- 
able publications to the Institute at the request of Harry F. 
Guggenheim. The Daniel Guggenheim [und ‘for the Promo- 
tion of Aeronautics made the publication of these digests pos- 
sible through a generous grant. They have been clipped and 
pasted on cards, and have made a valuable addition of ten 
thousand cards to the Index. 

The compiling of a card index of the Bibliography of the 
National Advisory Committee for Aeronautics is progressing 
and now includes all volumes from 1931—the latest published— 
to 1920. When this part of the Index is completed it should 
be possible to locate quickly all the papers and articles written 
by any author on one or more cards and to find all papers 
and articles written on any particular subject. It should be 
of the greatest value to anyone searching for aeronautical 
information. 

The Chronology of Aviation has been supplemented by chrono- 
logical information which is being received from many countries. 
Complete chronologies have been received from the Air Corps 
and the Weather Bureau. Other governmental agencies have 
such work in preparation. 

The Bibliographical file contains 8000 references which give 
information about the leading persons engaged in aeronautical 
activities. 

The Index will be as valuable to members of the Institute 
who do not live in New York as to those who can visit the 
Skyport. Inquiries regarding the books and articles on any 
subject will receive prompt attention. While the Index is still 
in a formative state it has been used by many persons with 
satisfactory results. 


many 


Les Fisches 


Lecture Tour 


Professor Albert Caquot, Director General of Services Tech- 
niques de l’Aeronautique of the French Air Ministry and 
Professor at L’Ecole National Superieure de 1’Aeronautique, 
will visit the United States on November 10th to begin a short 
lecture tour. He will speak, under the auspices of the Franco- 
American Institute of Science, on the subject “The mechanical 
characteristics of matter and contributions of French science, 
with recent applications to civil engineering and aviation.” 

Professor Caquot was educated at the Lycee of Reims and at 
the Polytechnical School of Paris. During the war he was 
Director of the Technical Section of Aviation and was an 
Officer of the Legion of Honor. 
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ENSTITUTE NOTES 


COMMENCEMENT OF THE INSTITUTE’S FourTH YEAR 


On October 1, 1935, the Institute will complete its third fiscal 
year. Reports as to the progress made during the year will be 
made by officers of the Institute at the Annual Meeting which is 
held in January. 

The membership of the Institute has grown during the year 
from 883 to 1141, an increase of 258 members. The financial 


condition has been sound owing to the strict economies prac- 
ticed in managing the affairs, particularly the expenses of the 
Journal which have been met by revenues from advertising. 
Over 97% of our membership have paid their dues for the 
year—an indication of the value placed on membership. 


A REQUEST 


The Institute is dependent on the prompt receipt of dues from 
members. Many accounts are payable at the beginning of our 
new fiscal year, October 1, 1935, when bills for annual dues 
are mailed to members. 

Members can show their appreciation of the valuable work 
done by the officers and staff in no more practical way than by 
paying their dues promptly on receipt of a bill. The Council 
has kept the dues lower than those of many other societies. 
Through contributions, the members have many facilities avail- 
able which would not be possible from the amount received 


from dues. 
Therefore, this request for immediate attention to this annual 


payment. 


PUBLICATIONS RECEIVED BY THE INSTITUTE 


The following aeronautical publications are received by the 
Institute : 


Aachener Abhandlungen 

Aero (Finland) 

Aero Digest 

Aeronautica Argentina 

L’Aeronautique (France) 

The Aeroplane (England) 

L’Aerotecnica (Italy) 

Les Ailes (France) 

Aircraft Engineering (England) 

Air Law Review 

Airways and Airports (England) 

Aviation 

L’Aviazione (Italy) 

Boletin de Aeronautica Civil (Argentina) 
Bulletin Technique du Bureau Veritas (France) 
Canadian Aviation 

Central Aero-Hydrodynamic Institute Reports (U.S.S.R.) 
Deutsche Luftwacht 

Envol (France) 

Les Fiches Aeronautiques (France) 

Flight (England) 

Indian Aviation 

Journal of Air Law 

Journal of the Institute of Engineers, Australia 
Journal of the Royal Aeronautical Society 
Journal of Scientific Instruments- (England) 
Journal of the Society of Automative Engineers 
Luftfahrtforschung (Germany ) 

Die Luftreise (Germany ) 

Model Airplane News 

Motor and Flying (Australia) 

National Aeronautic Magazine 

Official Aviation Guide 

The Pilot 

Popular Aviation 

Proceedings of the Institute of Radio Engineers 
Review of Physical Instruments 

Revista de Aeronautica (Spain) 

Revue du Ministere de I’Air (France) 

Rivista Aeronautica (Italy) 

Southwestern Aviation 

Sportsman Pilot 

U. S. Air Services 

Western Flying 


THE ENLARGED JOURNAL 


The Council, after receiving many requests to enlarge the 
page size of the Journal, decided that the reasons advanced for 
the change were sound, and commencing with the third volume 
have adopted what is known as the “standard” magazine size, 
The new page permits a use of larger drawings and illustrations 
and a greater elasticity in handling editorial matter. It also 
permits the Journal to use the same advertising plates which 
are prepared for other aeronautical publications. Heretofore, 
each advertisement for the Journal had to be especially prepared, 
photo engravings were multilated by reducing their size, and alf 
advertisements had to be specially handled by our printer at a 
heavy unnecessary expense to the Institute as well as much 
extra effort by friends of the Institute who support the Journal 
by advertising. 

Comments from members regarding any improvements which 
they believe will be helpful are always welcomed by the Editorial 
Board and the Council. 


Personnel Opportunities 


The Personnel Bureau serves organizations seeking to employ 
aeronautical specialists as well as individual members. The 
Bureau has been the means of arranging several very successful 
connections for members. 

Any member or organization may have their requirements 
listed without charge. 


Available 


Mechanical Engineer of wide experience as Executive and 
designer of power plants would connect up as Consultant or 
responsible head of design work or operation. Write Box 3l, 
Institute of the Aeronautical Sciences. 

Aircraft Engineer, sixteen years practical experience in all 
phases of work connected with complete design of Army, Navy 
and commercial airplanes. Capable of taking full charge of 
Engineering Department. Desires responsible connection with 
aircraft manufacturer. Would consider engineering connection 
with transport line or organization applying aircraft principles 
to non-aeronautical product. Write Box 32, Institute of the 
Aeronautical Sciences. 

College graduate, 27 years of age, with long training and 
experience along airline operations, traffic and maintenance 
work, and at present associated with an air line, desires a con 
nection as Field Manager or assistant to an air line executive 
Write Box 33, Institute of the Aeronautical Sciences. 

Junior Engineer, University graduate, with two years prac 
tical experience, one year stress analysis and aerodynamics, 
desires connection. Write Box 34, Institute of the Aeronautical 
Sciences. 

Junior Engineer, College graduate, with theoretical back- 
ground in Engineering Mechanics. Practical Experience if 
stress analysis and aerodynamics with government contractof 
Desires connection with engineering firm. Write Box 3§ 
Institute of the Aeronautical Sciences. 

Project and Design Engineer, 20 years experience on all 
phases of aeronautical engineering with leading American com 
panies. Capable of organizing and taking complete charge if 
either executive or engineering capacity. Also specialist if 
flight testing. Write Box 36, Institute of the Aeronautical 


Sciences. 
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